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NOMENCLATURE 


unit  vector  parallel  to  body  surface  in  plane  8  =  constant, 
positive  in  direction  of  increasing  x 

angle  from  meridian  0  =  constant  to  tangent  to  potential 
flow  streamline,  positive  in  direction  of  increasing  0 

constant.  Equation  (40) 

qe  -  ur0sin  a 

element  of  area  on  resultant  body,  body  plus  displacement 
surface 

c  x  a,  unit  vector  parallel  to  body  surface  in  plane 
x  =  constant 

constant,  Equation  (40) 

u>r0sin  a 

unit  vector  normal  to  body  surface,  outward  from  surface 
constant,  Equation  (79)  *  / 

constant.  Equation  (80) 


c2/Ci 

constant.  Equation  (83) 
tur0cos  a 

drag  coefficient  based  on  area  nrM2 

Magnus  yawing  moment  coefficient  (Equation  (154)) 

normal  force  coefficient  (Equation  (152)) 

P-p® 


pressure  coefficient 


Poo  ~  2 
2~  V°° 


Magnus  side  force  coefficient  (Equation  (153)) 

displacement  of  origin  of  x,  r,  0  system  from  origin  of 
£,  h,  8  system  (Figure  8) 


vi 


/ 


NOLTR  72-80 


unit  vector  .  direction  of  free  stream  velocity  (Figure  8) 
unit  vector  normal  to  and  in  x  plane  (Figure  8) 

£>  x  ^2  (Figure  8) 

Fy  +  i  Fz,  complex  Magnus  force  on  portion  of  body  between 
crors  section  %  =  C  and  b  jdy  nose 

)mponeftw  of  Magnus  fo^ce  F  along  Y  axis  (Figure  6) 

■  Dmp*'rarv  •>/  Magnus  corce  F  along  L  axis  (Figure  6) 

^ anctic  .  " r  n  (See  Table  I) 

radial  distance  in  plane  £  =  constant  (Figure  8) 

rau; il  distance  to  point  on  body  in  plane  £  =  constant 
(Figure  6) 

unit  radial  vector  in  plane  £  =  constant  (Figure  8) 
unit  vector  along  body  axis  (Figure  8) 
vd 

radius  of  starting  circle,  (Equation  (100)) 

A 

unit  vector  normal  to  i  and  in  x,  £  plane  (Figure  8) 
parameter  (See  Table  III) 

.constant  in  friction  formula  (Equation  (95)) 
parameters  for  ellipsoid  (Equations  (150)  and  (151)) 
exponent  (Equation  (40)) 
reference  length 

exponent  in  friction  formula  (Equation  (95) ) 

Mach  number 

Magnus  moment  about  nose  of  body  acting  on  portion  of 
body  between  cross  section  £  =  C  and  the  body  nose, 
positive  as  shown  in  Figure  1 
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exponent  in  velocity  profile  formula  (Equation  (52)) 

unit  vector  normal  to  body  plus  displacement  surface, 
outward  from  surface 

number  of  base  radii  behind  base  beyond  which  wake 
thickness  is  constant  and  equal  to 

vector  normal  to  body  plus  displacement  surface,  outward 
from  surface 

— P" 

spin  rate  parameter  — 

Vco 

static  pressure 

velocity  at  outer  edge  of  boundary  layer 
right  hand  side  of  Equation  (24) 
velocity  vector  (Equation  0.05} ) 
radial  distance  in  plane  x  =  C  (Figure  8) 
radius  of  body  of  revolution 
maximum  radius  of  body  of  revolution 

radial  coordinate  of  section  of  resultant  body  in  Z'  plane 
(Figure  7) 

unit  vector  along  radial  direction  in  x  =  constant  plane 
(Figure  8) 

radial  coordinate  of  resultant  body,  rQ  +  6* 

reference  Reynolds  number  V°°L 

v 

displacement  radius  of  wake  (Equation  (C-l) ) 

constant  displacement  radius  of  wake  far  behind  body 
(Equation  (C-l)) 

distance  on  bods  surface  along  potential  flow  streamline 
distance  along  body  surface  in  plane  8  =  constant 
distance  along  body  surface  in  plane  x  “  constant 

distance  along  section  of  body  plus  displacement  surface 
in  plane  £  =  constant  (Figure  7) 
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vector  from  origin  of  £,  h,  6  system  of  coordinates 
(Figure  8) 

thickness  ratio,  minor  axis  divided  by  major  axis  for 
ellipsoid 

temperature 

velocity  parallel  to  surface  and  in  plane  0  =  constant 

J.. 

value  of  Q  e  at  stagnation  point 
dsx 

velocity  parallel  to  surface  and  in  direction  of  q& 

velocity  parallel  to  surface  and  in  plane  x  =  constant 

velocity  parallel  to  surface  and  normal  to  direction  of 

% 

velocity  along  circle  in  Z  plane  (Figure  7) 

velocity  along  cross  section  in  Z'  and  zq  planes  (Figure  7) 

velocity  normal  to  circle  in  Z  plane  (Figure  7) 

velocity  normal  to  cross  section  in  Z '  and  zQ  planes 
(Figure  7) 

magnitude  of  free  stream  velocity 

velocity  vector  along  displacement  surface  (Equation  (10)) 
velocity  normal  to  body  surface 
distance  along  axis  of  revolution 

value  of  x  at  intersection  of  axis  of  revolution  and 
plane  K  =  C 

distance  from  body  surface  in  direction  normal  to  surface 
value  of  y  slightly  larger  than  6 
coordinate  in  plane  (Figure  6) 

complex  plane  coinciding  with  plane  £  =  constant  and 
having  origin  of  coordinates  on  £  axis,  zQ  =  Y  +  i  Z 
(Figure  6,  7) 

complex  plane  coinciding  with  plane  £  =  constant  and 
having  origin  of  coordinates  on  x  axis  (Figure  7) 


NOLTR  72-80 


Z  complex  plane  in  which  section  of  body  plus  displacement 

surface  is  a  circle  as  a  result  of  a  transformation 
(Figure  7) 

a  angle  of  attack  of  body 

0O  tan-1  i p 

$  angular  coordinate  in  ^  =  constant  plane  (Figures  6,  7,  8) 

A 

0  unit  vector  in  direction  of  increasing  0  (Figures  6,  8) 

Y  angular  coordinate  in  Z'  plane  (Figure  7) 

A  integration  interval  2 tt  divided  into  equal  steps  of 

length  A  (Equation  (144)) 

6  boundary  layer  thickness 

6*  boundary  layer  displacement  thickness  defined  by  Equation 

(23) 

7  gradient  symbol 

£  (o  -  y)  see  Equation  (128) 

c  y/s 

n  starting  angle  (Equation  (100)) 

0  unit  vector  normal  to  r  and  in  direction  of  increasing 

9  (Figure  8)  (0  =  r  x  I) 

0  angular  coordinate  (Figure  2) 

k  ratio  of  specific  heat  at  constant  pressure  to  specific 

heat  at  constant  volume 

X  radial  coordinate  in  circle  plane  (Figure  7) 

p  viscosity 

v  kinematic  viscosity  u/"p 

A 

v  unit  vector  normal  to  section  of  body  £  =  constant 

(Figures  6,  7) 

£  distance  in  direction  of  free  stream  velocity  (Figures  6,  8) 

p  density 

a  angle  in  Z  plane  (Figure  7) 
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w 


Twx 


lw. 


wall  shear  stress 

wall  shear  stress  along  0  =  constant  direction 

wall  shear  stress  along  x  =  constant  direction 

angle  between  surface  of  body  and  x  axis  in  plane 
9  =  constant  (Figure  2) 

three-dimensional  perturbation  velocity  potential 
(Equation  (105)) 

two-dimensional  perturbation  velocity  potential  in  plane 
5  =  c 


angle  defined  by  Equations  (12)  and  (13) 


*2 


angle  defined  by  Equations  (14)  and  (15) 


ip  tangent  of  angle  measured  from  direction  of  potential 

flow  streamline  to  direction  of  surface  shear  stress, 
positive  in  direction  of  increasing  0 

a)  angular  spin  velocity  of  body  -r-—*,  positive  as  shown  in 

(Figure  1)  v<» 

fi  defined  by  r'  =  rQefi  (Equation  (124)) 

Qq  defined  by  A  =  r0e^o  (Equation  (125)) 


Subscripts 

b  at  base  of  body 

e  at  outer  edge  of  boundary  layer 

M  maximum  value 

s  at  stagnation  point 

w  at  surface 

W  wake 

6*  on  displacement  thickness  surface 

very  far  ahead  of  body 
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Superscripts 

!!  fluctuating  quantity 

<  >  mean  value 

*  value  very  far  behind  base 

_  dimensional  quantity 
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INTRODUCTION 


When  a  spinning  body  of  revolution  is  flying  at  an  angle  of 
attack,  a  force  normal  to  the  angle  of  attack  plane  acts  on  the  body. 
Associated  with  this  side  force  is  a  moment.  The  force  and  moment 
are  known  as  the  Magnus  force  and  moment.  Although  the  Magnus  force 
is  usually  only  a  small  fraction  of  the  normal  force  it,  and  its 
moment  about  the  body  center  of  gravity,  can  have  an  important  effect 
on  the  body  trajectory.  In  order  to  predict  the  Magnus  force  and 
moment  and  in  order  to  better  understand  experimental  results  a 
theory  is  needed. 

The  problem  is  to  calculate  the  Magnus  force  and  moment  given 
only  the  shape  of  a  body  of  revolution,  its  speed  along  the  flight 
path,  its  spin  rate,  its  angle  of  attack,  and  the  properties  of  the 
atmosphere.  If  the  flow  were  symmetric  about  the  plane  formed  by 
the  axis  of  the  body  and  the  free-stream  velocity  vector  there  would 
be  no  side  force,  a  force  perpendicular  to  the  plane  of  symmetry. 
There  would  also  be  no  yawing  moment.  Because,  however,  the  fluid 
through  which  the  body  is  moving  is  viscous,  a  boundary  layer  is 
present  on  the  body.  The  spin  of  the  body  combined  with  the  angle 
of  attack  causes  the  boundary  layer  to  be  unsymmetric  with  respect 
to  the  plane  of  symmetry,  the  plane  formed  by  the  axis  of  revolution 
and  the  free-stream  velocity  vector.  The  resultant  configuration 
is  shown  in  Figure  1.  Because  the  symmetry  of  the  flow  about  the 
angle  of  attack  plane  is  destroyed  by  the  unsymmetric  boundary  layer 
a  side  force  and  moment  exist.  For  the  symmetry  to  be  destroyed, 
both  spin  and  angle  of  attack  must  be  present. 

In  the  present  investigation  the  Magnus  force  and  moment  are 
calculated  for  a  body  with  no  separation  of  the  boundary  layer. 

The  boundary  layer  is  therefore  thin  over  the  entire  body  and  there 
are  no  regions  of  vorticity  shed  into  the  flow  ahead  of  the  body 
base.  To  calculate  the  Magnus  force  and  moment  for  such  a  flow  the 
inviscid  flow  around  the  body  is  found  first.  Then  the  boundary 
layer  displacement  thickness  surface  surrounding  the  body  is  cal¬ 
culated  by  use  of  boundary  layer  theory.  The  displacement  thickness 
surface  is  added  to  the  body  and  the  force  and  moment  calculated  for 
the  resultant  body  in  an  inviscid  flow. 

This  method  for  the  calculation  of  the  Magnus  force  and  moment 
is  that  of  Martin  (Reference  1)  who  presents  a  theory  for  the  Magnus 
force  and  moment  on  a  cylinder  at  a  small  angle  of  attack.  The 
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boundary  layer  is  laminar  and  incompressible.  Quantities  of  higher 
order  than  the  second  in  angle  of  attack,  spin  velocity,  and  distance 
from  the  leading  edge  and  their  products  with  one  another  are  neglect¬ 
ed.  Martin  finds  that  the  Magnus  force  is  directly  proportional  to 
the  product  of  the  angle  of  attack,  spin  velocity  ratio,  body  length, 
and  the  displacement  thickness  at  the  body  base  at  zero  angle  of 
attack.  For  turbulent  flow  Martin  replaces  a  calculated  constant 
of  proportionality  by  an  unknown  coefficient. 

Platou  (Reference  2)  presents  experimental  data  for  supersonic 
flow  and  concludes  that  Martin's  incompressible  flow  theory  predicts 
the  correct  order  of  magnitude  for  3  to  5  caliber  bodies  with  laminar 
boundary  layers.  To  apply  Martin's  theory  to  a  bullet  shaped  body 
an  arbitrary  allowance  must  be  made  for  the  nose  portion  of  non¬ 
constant  diameter.  The  extrapolation  of  Martin's  theory  to  turbulent 
flow  correctly  predicts  the  Magnus  force  to  be  directly  proportional 
to  the  product  of  angle  of  attack  and  spin  rate.  Platou  also  finds 
that  the  theory  of  Kelly  and  Thacker  (Reference  3),  which  includes 
a  radial  pressure  gradient  and  skin  friction  effect,  does  not  agree 
with  Martin's  prediction  nor  does  it  agree  with  the  experimental 
result  that  the  Magnus  force  depends  on  the  spin  to  the  first  power. 

Sedney  (Reference  4)  calculates  the  Magnus  force  and  moment  on 
a  slender  spinning  cone  at  a  small  angle  of  attack  in  supersonic 
flow.  The  boundary  layer  flow  is  laminar.  Terms  of  higher  order 
than  the  first  in  angle  of  attack,  spin  velocity,  and  in  the  product 
of  spin  velocity  and  angle  of  attack  are  neglected.  The  method 
follows  Martin  and  also,  like  Martin,  Sedney  uses  slender-body  theory 
to  calculate  the  force  and  moment  on  the  body  that  results  when  the 
boundary  layer  displacement  surface  is  added  to  the  body  of  revolution. 
No  experimental  test  of  the  predictions  seems  to  be  available. 

In  the  present  investxgation  a  method  is  developed  to  calculate 
the  Magnus  force  and  moment  for  a  body  of  revolution  of  general 
shape  with  an  unseparated  turbulent  boundary  layer.  Although  the 
boundary  layer  calculation  is  for  incompressible  flow,  the  calculation 
method  can  be  used  for  Mach  numbers  up  to  the  transonic  range  because 
Mach  number  effects  on  boundary  layer  flow  are  usually  small  for 
local  Mach  numbers  less  than  unity.  The  boundary  layer  calculation 
is  based  on  the  momentum  integral  method.  The  force  on  the  body  is 
calculated  by  slender-body  theory  (References  5  and  6).  A  discussion 
is  given  of  the  calculation  of  the  force  by  a  more  exact  method 
than  slender-body  theory,  namely,  the  method  of  Hess  and  Smith 
(References  7  and  8) . 


ANAYLSIS 


In  order  to  calculate  the  force  and  moment  on  a  spinning  body 
the  effective  shape  of  the  body  in  an  inviscid  flow  is  needed.  The 
effective  shape  is  found  by  calculating  the  boundary  layer  displace¬ 
ment  surface  and  adding  it  to  the  body  of  revolution.  The  displace¬ 
ment  surface  6*  is  the  surface,  which,  added  to  a  body  in  an  inviscid 
flow,  results  in  a  body  with  the  same  streamlines  as  those  outside 


2 


NOLTR  72-80 


the  boundary  layer  in  the  real  viscous  flow  (Reference  9).  Therefore, 
because  according  to  boundary  layer  theory  the  pressure  difference 
across  the  boundary  layer  is  negligible,  the  pressure  on  the  real 
body  in  the  viscous  flow  is  the  same  as  on  the  body  plus  5*  in  the 
inviscid  flow. 

Displacement  Surface 

To  find  the  displacement  surface,  6*,  use  is  made  of  the  fact 
that  the  streamlines  outside  the  boundary  layer  are  not  changed  by 
replacing  the  body  with  the  boundary  layer  over  it  in  the  viscous 
flow  by  the  body  with  6*  added  to  it  in  an  inviscid  flow.  Conse¬ 
quently  the  velocity  component  w  on  these  streamlines  is  also  unchang¬ 
ed.  To  find  w  on  a  streamline  near  the  outer  edge  of  the  boundary 
layer  the  continuity  equation  is  used;  for  a  steady  flow  with  the 
coordinate  system  shown  in  Figure  2  it  is 


cos<f> 


3__ 

9x 


(pur)  +  f^-(pv)  +  fy  (pwr)  =  0 


(1) 


where  p,  u,  and  v  are  time  mean  values  and  pw  is  (pw  +  <p"w">). 
Equation  (1)  is  the  same  as  Equation  (19)  page  414  Reference  10  when 
r  is  put  equal  to  rQ.  The  coordinate  system  is  fixed  in  the  fluid 
and  the  body  rotates  around  its  axis  of  symmetry.  All  quantities 
are  non-dimensional;  the  velocities  are  non-dimens ionalized  by 
the  lengths  by  IT,  the  densities  by  p^,  and  the  pressures  and  shear 
stresses  by  p  00  Voo  . 

An  integration  of  (1)  with  respect  to  y  up  to  yq,  where  yq  is 
slightly  greater  than  the  boundary  layer  thickness,  results  in 


|j(pur)dy  -  | 


36 


(pv)dy 


J 


o 


J 


o 


(2) 
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For  an  inviscid  flow  over  the  body  with  6*  added  to  it  the  result  for 
( pW )  y q  is 


c  os4> 
r 


cix (Peuer) dy  -  r  fe(peVdy  +  {pw)6* 


The  term  (pw)  appears  because  the  6*  surface  has  a  slope  with 
respect  to  the  body  and  w  is  normal  to  the  surface  without  6*. 

In  Equation  (3)  p,  u,  and  v  have  the  subscript  "e"  because  the  flow 
is  inviscid  and  the  velocity  varies  negligibly  slowly  with  y. 

Then,  because  the  streamline  at  y  is  the  same  for  both  the 
viscous  and  the  inviscid  flow  pw  is  also  the  same.  Therefore, 
equating  (2)  and  (3)  the  result  is 


COS<P 

r 


g^[r(peue-pu)]dy  -  -  -(p^-pvjdy 


(P  u  r)dy  +  — 
9x  e  e  J  r 


ae(ceve)dy  +  (ou)6* 
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-  cosft  3_ 


r  3xo 
o 


(Peue-pu)dy  -  !j 


r 


J 


(P  eV  e"pv)dy  +  F"cos<j>fx(ropeUe) 
e  e  o 


+  r”  fe(peve>  +  (pw)6*  =  0 
o 


(9) 


To  find  (pw)^*  use  is  made  of  the  fact  that  in  the  inviscid 
flow  over  the  body  with  6*  added  to  it  the  velocity  component^normal 
to  the  displacement  surface  6*  is  2ero.  The  velocity  vector  Ve 
along  the  6*  surface  is 


V  =  a  u  +  b  v  +  cw  (10) 

e  e  e  e 


where  a  is  a  unit  vector  along  the  body  surface  and  in  the  x 
direction,  B  is  a  unit  vector  along  the  body  surface  in  a  plane 
x  =  constant,  and  Q  is  a  unit  vector  normal  to  the  body  surface  and 
outward.  The  condition  for  no  velocity  normal  to  the  displacement 
surface  is 


V 

e 


N  = 


0 


(11) 


where  N  is  a  vector  normal  to  the  body  plus  6*  and  outward  from  the 
surface.  To  find  N  note  that  ft  is  normal  to  each  of  the  two  vectors 
(a  cos  Xi  +  c  sin  Xi>  and  (6  cos  X2  +  c  sin  X2)  which  lie  in  the 
displacement  surface  along  0  =  c  and  x  =  c, respectively .  The  angles 
Xi  and  X2  are  gi^en  by 


cosxl  ■  5^  (12) 
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sinxL 


As 

3gK.i!* 


cosx2 


r  A0 
o 


sinx2 


36* 

30 


A0 


where 


D,  = 


T 


As„  +  ( 


x 


^As  >2 

9sxisx> 


and 


D2  ^ 


r  2A02 
o 


(|fA0)2 


Then 


N 


{a  cosx-l  + 


-  sinx^)  x  (b  cosx2  + 


c  sinx2) 
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Equation  (23)  is  the  equation  for  the  displacement  surface  6*  (see 
Reference  9) .  The  Equation  (23)  is  a  partial  differential  equation 
for  6*.  When  Equation  (23)  is  expanded  the  result  is 


36*  .  ,  36* 

r  cosAp  u  t- —  +  p  v  tt-q— 
o  e  e3x  e  e30 


-  6*  [rocos4>  -3x 


3peue  9ro 

+  cos^-  peue  + 


dpeve. 


+  cos<|>3~  f  <peVpu)dy  +  rocos<|lfx  J  (peU 


-pu) dy 


vj 

I of  (peV 


pv)  dy 


(24) 


where  the  upper  limit  yq  in  (23)  has  been  replaced  by  6  in  (24); 

6  is  the  smallest  value4of  y  for  which  simultaneously  pu  =  peuQ  and 
pv=peve. 

Equation  (24)  is  a  Lagrange  linear  partial  differential  equation. 
To  obtain  a  solution  the  two  subsidiary  equations 


dx 


d0 


r  cos<j>p  u 
o  e  e 


(25) 


P  v 
Fe  e 


and 


dx 


d6  5 


rocos4»peue 


(26) 
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are  to  be  integrated;  the  quantity  Q  is  the  right  hand  side  of 
Equation  (24).  Equation  (25)  is 

♦ 


de 

dx 


v 


e 

rQUecos<}> 


(27) 


and  is  the  equation  of  an  inviscid  flow  streamline  over  the  body.  To 
integrate  (26)  requires  that  the  integrals  in  (24)  be  known.  That  is, 
the  boundary  layer  velocity  distribution  and  thickness  must  be  known. 
When  the  boundary  layer  properties  are  known,  Equations  (25)  and 
(26)  are  integrated  together  thereby  giving  6*  along  a  streamline. 

By  integrating  (25)  and  (26)  along  a  sufficient  number  of  streamlines 
the  displacement  surface  6*  on  the  body  is  calculated. 

Momentum  Integral  Equations 

For  thin  boundary  layers  the  equations  of  motion  are  (Reference  10 
page  14) 


,3u  .  pv  3u  .  3u  cos<|> 
Pucos^  +  ~  Tfl  +  Pw  ~  T- 


r  30 
o 


3y 


3r 

i 

3>c 


o  2 

pv  = 


,3p  9Tx 
'  COS<^3x  +  37" 


(28) 


and 


.  3  v  , 
pucos<J>—  + 

dX 


pv  3v 

r  30 
o 


,  3r 

cosy  _ o 

r  3x 
o 


puv 


1  3P  3t0 
r  36  +  ay 


(29) 


and  the  continuity  equation  is 


cos4>  3  ,  .13  a 

r_  3x(puro}  +  FJ  31T(pv)  +  3^(pw)  =  0 


(la) 
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To  obtain  the  momentum  integral  equations  the  x  equation  of  motion 
(28)  and  the  0  equation  of  motion  (29)  are  integrated  with  respect 
to  y  in  the  usual  way  with  the  help  of  the  equation  of  continuity. 
The  result  for  the  x  equation  is 


roc°s^ 


6  <$  ^ 

5  I  pu^-u.dy  +  sin*  /  pu(ue-u)dy  +  rocos<j,^  J  (p^ 


(30) 


-pu)dy 


a 

g  f  9  U 

+  36  I  PvtUg-uldy  *  JT  I 


j 


-pv)dy  -  sin <J>  /  (p  v  2-pv2) dy  =  t  r 
'  e  e  J  wo 

x 


For  the  0  equation  the  result  is 


3v 

IF  /  pv(ve-v)dy  +  35-  /  (Peve-pv)dy  *  r0cos«,£ 


/ 


/ 


x  I  -(Vv)dy 


( 

/ 


9v 


+  sin<}:  /  pu(v  -v)dy  +  r  cos<j>  ^ —  /  (p  u  -pu)dy  +  sind>  /  (p  u  v  -puv)dy 
I  e  o  :X  I  o  &  /  e  e  e 


=  t  r 
w  >  o 


(31) 


These  equations  are  the  same  as  on  page  416  of  Reference  10  in 
different  notation.  In  order  to  use  the  momentum  Equations  (30) 
and  (31)  to  calculate  the  terms^/1'  (peue-pu)dy  and q/6  (peve~pv)dy 

in  Equation  (24)  for  the  displacement  surface  5*  an  expression  for 
the  velocity  profile  through  the  boundary  layer  is  needed. 
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Velocity  Profile 

The  velocity  profile  of  the  velocity  component  in  the  direction 
of  the  outer  streamline  is  assumed  to  be 


Y7Jfrw  „  f  (?)  (32) 

qe  q,w 


where  ?=y/6 ,  and  the  velocity  profile  for  the  velocity  component 
perpendicular  to  the  outer  streamline  is  assumed  to  be 


■SL_-..a.*.S£  =  y g  (C)  f (C)  (33) 

■e  q,w 


In  (32)  and  (33)  qe  is  the  speed  along  the  streamline  at  the  outer 
edge  of  the  boundary  layer;  uq  is  the  speed  in  the  direction  of 
the  outer  streamline  and  vq  is  the  speed  ir.  the  direction  perpendic¬ 
ular  to  the  outer  streamline.  The  quantity  ip  is  an  as  yet  unde¬ 
termined  parameter.  At  the  surface  (see  Figure  3) 


'v.  ,,  =  wr  sina, 
q  ,w  o 


(34) 


and 


v  =  wr  cosa. 
q,w  o  / 


(35) 


where  a  is  the  angle  between  the  outer  streamline  and  a  line  0=c. 
Also,  from  (32)  f(o)=0  and  f(l)=l.  At  ?=1(  vq=0  by  definition. 
Therefore  from  (33)  and  f  ( 1 ) =1  it  follows  that 


9  (1) 


-  v 


rW 


(36) 
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Also,  division  of  (33)  by  (32)  results  in 


Sj. 


w 


u 


q,w 


vg 


(37) 


As  y->-o  (37)  becomes 


1 

y-*o 


u 


-  u 


q,w 


g  (o) 


(38) 


The  quantity  g(o)  is  put  equal  to  unity.  Then  (38)  becomes 


1 

y-*o 


(39) 


3  v 

_  q 

The  quantity  y  (3 — )  equals  the  shear  on  the  surface  in  a  direction 


r4  \ 

w  3- 


w 


normal  to  the  outer  streamline;  a  positive  value  means  the  shear  acts 

3Hq 

in  the  direction  of  increasing  8.  The  quantity  yw( - )  is  equal  to 

^y  w 

the  shear  acting  on  the  surface  in  the  direction  of  the  outer  stream¬ 
line  flow.  Consequently  ’•>  is  the  tangent  of  the  angle  measured  from 
the  direction  of  the  outer  streamline  to  the  direction  of  the  surface 
shear  stress;  ^  is  positive  in  the  direction  of  increasing  6. 


The  function  g (?)  is  taken  to  be 


gU)  =  U-C)£  +  ao  +  bQ  e 


(40) 
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The  condition  g(o)=l,  previously  imposed,  gives  ao=0..  From  the 
expression  (36)  for  g(l)  and  from  (40),  the  result  is 


b 


o 


♦'Ww1 


Therefore  (40)  becomes 


g(?) 


(1-0  - 


. 

♦“Ww1 


(41) 


The  exponent  l  is  found  by  considering  conditions  at  the  stagnation 
point  of  a  rotating  body  at  an  angle  of  attack ;  thus  the  stagnation 
point  is  not  on  the  axis  of  rotation.  At  the  stagnation  point  the 
outer  stream  velocity  is  zero  and  the  surface  has  a  rotational 
velocity  wr0.  Let 


v  =  wrQF(?) 


(42) 


where  F(o)=l  and  F(1)=0,  that  is,  the  fluid  has  velocity  tor0  at  the 
surface  and  zero  velocity  outside  the  boundary  layer.  Also,  u,  the 
velocity  in  the  direction  of  sx  is  zero  at  the  stagnation  point. 

The  relations  between  the  velocity  components  u_  and  vq  and  the 
velocity  components  u  and  v  are 


u  =  ucosa  +  vsina  (43) 

q 


and 


v  =  -  usina  +  vcosa  (44) 

q 
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At  the  stagnation  point  (43)  and  (44)  becomes,  with  (42), 

Ug  =  uiXQ  F  sina  (45) 

and 

v  =  wr  F  cosa  (46) 

Hi 

Equations  (45),  (34),  (32)  and  qe=0  result  in 

F  =  1-f  (47) 

From  (34),  (35),  qe=0,  (46)  and  (47)  are  used  in  (33)  the  result  is 

cota  =  i}/  g  (48) 

From  (39),  (45)  and  (46)  it  follows  that 

cota  =  ip  (49) 

Therefore  from  (48)  it  follows  that  g=l  for  all  c  at  the  stagnation 
point.  Now  consider  Equation  (41).  At  the  stagnation  point  qe=0, 
ip-cot  a,  and  relations  (45)  and  (46)  with  F(o)=l  hold. 

Equation  (41)  then  becomes 

gU)  =  (1-0*  +  5  (50) 
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Therefore  to  have  g(C)~l  at  the  stagnation  point,  put  £=1. 
Then  (41)  becomes 


The  velocity  profile  (32)  is  taken  to  be  a  power  profile  by 
putting 


f(C)  =  ?n 


(52) 


Because  the  velocities  tnat  appear  in  the  two  momentum  integral 
Equations  (30)  and  (31)  contain  u  and  v  it  is  necessary  to  express 
u  and  v  in  terms  of  the  velocities  Uq  and  Vq  that  appear  in  (32) 
and  (33) .  The  relations  are 


u  =  u  cosa-v  sina  (53) 

q  q 


and 


v  =  u  sina  +  v  cosa  (54) 

q  q 


When  Uq  and  Vq  are  found  from  (32)  and  (33)  and  (34) ,  (35),  (51), 
and  (5?)  are  used.  Equations  (53)  and  (54)  become 


u  =  (A£n+B)  cosa  -  [A^(Cn-Cn+1)  +  C(l-Cn+1)]  sina  (55) 
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and 


v  =  (Asn+B)sina  +  [A4’(Cn-?n+1)  +  C(l-Cn+1)]  cosa 


(56) 


where 


A  = 

B  = 

C  = 


a  -u  =  q  -tor  sina 
4e  q,w  ^e  o 


u  =  tor  sina 
q,w  o 


v  =  tor  cosa 
q,w  o 


(57) 

(58) 

(59) 


Differential  Equations  for  6  and  t|j 

Expressions  (55)  and  (56)  are  put  into  the  integral  momentum 
Equations  (30)  and  (31)  and  the  integrals  evaluated.  The  density  is 
taken  constant  from  here  on  in  the  boundary  layer  analysis.  A 
sample  term  is  _^(ue-u)dy  which  is  written  as  5 (ue-u) 

The  result  is 


] 

/ 


„  1.  ] 

(ue-u)d£  =  A  ^—j-cosa  -  A^j~jsina  +  (A^+O^-sina 


or 


(ue-u)d^=  AG^cosa  +  !|;AG2sina  +  CG^sina 
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where 


c  _  1 
2  (n+1)  (n+2) 

r  -  n+1 
G4  "  n+2 


-  n 
n+1 


There  are  fourteen  G's;  they  are  listed  in  Table  1. 

In  Table  II  are  listed  the  integrals  that  occur  in  Equations  (30) 
and  (31) .  The  procedure  results  in  two  partial  differential  equations 
with  6  and  ip  as  dependent  variables  and  x  and  0  as  independent  vari¬ 
ables.  The  x  equation,  the  result  of  (30)  and  (55)  and  (56)  is 


rQcos4,|^(l)i2j1+-|iJ2+J3J  +  rQCos$|i  «[2*J1+J2] 


+  Hi*2j6+,|’J7+J81  +  H  M2+J6+J7J 


=  ro  Twv  *  S(I*  S3T  +  *TT  +  3riro“s* 


3r 


3u 


+  [1»  J1+i('J2+J3]g^-cos<})  +  roCOS({)y^-[lj;J4+J5] 


^3»J  £  3  3*Jp  3u 

+  [^  dir*  3F"J  +  30“^J18+J19^ 


(60) 


-  3r 

-  [^^J9+^J10+J11]^cos4.} 
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The  8  equation,  the  result  of  (31)  and  (55)  and  (56)  is 


roCOs4,H[,|;2j6+^J14+J15]  +  roCOS'0  [2^VJ14] 


+  H^2j9+1^J12+J13]  +  II  6t2^J9+Ji2] 


(61) 


_  3J,.  «  J-,  ,  8J-i  c  2  ^ro 

-  v,(  -  +  3F^roC°S«>  +  t*  J6+*J14+J15113TCOS't' 


3v 

+  g~  [(|/J4+J5]rocos<|)  +  [t|> 


2  3J9 


38  +  ^36 


3J12  ,  3J131 


38 


Sv  2  *r_ 

r  ..  -r  .  -*•  1  .  f  .1  .*•  • 


+  JT  [*J18+J19]  +  [*  J6+^ J16+J17 1  Sx2®08* ] 


The  J's  are  listed  in  Table  III.  They  are  independent  of  6  and  ip 
and  are  calculated  from  the  inviscid  flow  velocity  distribution  over 
the  body,  spin  rate,  and  the  distribution  rD(x). 

The  simultaneous  solution  of  (60)  and  (61)  gives  |A  AA,  AA 

3x  98  3x 


and  || 

and  (61)  are  solved  for  (|A)  and  (|A)  with  |A  =  o  and  1A  =0.  Then 


jg  at  a  point  x,0.  An  iteration  process  is  used.  First  (60) 


1  3x  x  38  ~  30 

the  values  of  (AA)  and  (|A)  are  substituted  into  (60)  and  (61),  and 

dX  dX  ^ 

(60)  and  (61)  are  solved  for  (|A)  and  (||)  .  Then  (iA)  and  (iA) 

u  1  oi)  j  382. 

are  substituted  into  (60)  and  (61),  and  (60)  and  (61)  are  solved  for 

(^A)  and  (|A)  .  These  values  of  (^A)  and  (Aii)  are  substituted  into 
<?x  2  ^x  2  ox  2  3x  2 

(60)  and  (61)/  and  (60)  and  (61)  solved  for  (~)  and  (^-)  ♦  These 

O  0  a  U  n 

'NT  „  Z  Z 

values  are  then  used  to  find  ill )  and  (lx.)  which  are  used  to  find 

oX  3  :-x  3 

(11.)  and  (®V)  and  so  on  until  further  iteration  produces  a  neglig- 


‘30 


l?0' 


3  ou  2 

ible  change  in  (^)  t  (11)  t  (11)  and  (14)  .  in  the  calculations  the 
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-  1 


iteration  is  stopped  either  after  25  iterations  or  when  t-* -  . 

^  si-l  ;x+Ax 

and  (— —  -1)  are  both  less  than  .001.  The  subscript  i  is  the 

^i-1  x+Ax 

iteration  number;  6x+^x  and  ’^x+^x  are  the  values  of  6  and  \p  at  the 
next  point  along  a  streamline.  They  are  found  from 


A<5 


+  3  0  dx  1  AX 


0 


(62) 


and 


**  ’  fe>  ♦  (H>  I 

i*  U  x 


Ax 


(63) 


where  is  obtained  from  Equation  (27),  the  equation  for  a  streamline, 
dx 

Initial  Conditions  for  6  and  ip 

In  order  to  begin  the  integration  of  (60)  and  (61) ,  the  values 
of  6  and  at  the  stagnation  point  are  needed.  Up  to  now  a  method 
for  calculating  6  at  the  stagnation  point  of  a  spinning  body  when 
the  stagnation  point  is  not  on  the  axis  of  rotation  has  not  been  found. 
Consequently,  the  needed  value  of  6  is  obtained  by  extrapolation  from 
the  value  of  S  at  the  stagnation  point  on  a  spinning  body  at  zero 
angle  of  attack  and  the  value  of  <5  at  the  stagnation  point  of  a  non¬ 
spinning  body  at  an  angle  of  attack.  Thus,  it  is  assumed  that  6  at 
the  stagnation  point  when  spin  and  angle  of  attack  are  both  present 
can  be  written  as 


6 


a  ,p 


6  + 

0,0 


(!£) 

^3p; 


P  + 


0,0 


i2i) 


0,0 


(64) 
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At  a  =  0  the  relations  u  =  u'  r  ,  and  r  =  s  hold.  Also 

e  e ,  s  o  o  x 


inr1  oos  ♦  = 


The  friction  coefficient  is  written  as 


—  ,3ux 


-  3u/ue 

*  hr?—: 1 


U  2  "p  u  2  V  u  6L 

e  e  00  e 


9u/u 

/ - 

1  9?  \ 


r  6  Re_ 
e/S  o  L 


When  these  relations  are  substituted  into  (68)  ,  the  result  is 


96G?Gg 


u'  r  6 
e,s  o 


9u/u 


/  2 
-  u  _6^ 


e  ,s  [3G7G8+G7-^^Gi;l 


+  ats_  q  _ 


The  quantity  SG^Gg  is  equal  to  the  ratio  of  the  momentum  thickness  to 
the  boundary  layer  thickness  6,  fl  Cn(l-£n)dC.  At  the  stagnation 


point  rD  =  0.  Consequently  to  avoid  an  infinite  value  of 


9  SGyGg 
3s  ' 
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which  is  physically  unrealistic,  it  is  necessary  that  the  term 
inside  the  braces  on  the  right  hand  side  of  (70)  be  zero;  that  is. 


3u/u 

( _ _ e,) 

'  as  * w 

,2  /  „ 

6  u  ^ReT 
©  z  s  L 


(3G?G8+G7)  +  trGn  - 


u 7  5  T5- 

e,s  e,s 


When  the  relations  (66)  and  (67)  are  used  in  (61)  and  the  same 
procedure  used  that  was  used  to  obtain  (71)  the  result  is 


9v/ue 

The  ratio  ( — *-— )  can  be  replaced  by  an  expression  involving 

3u/ue  3  5  w 

( - )  by  using  (56)  with  "a"  =  0.  That  is 

3?  w 


v  =  ueiMCn-?n+1)  +  u>rQ(l-Cn+1)  (73) 


Then 


— =  ij>  [n^11**1-  (n+1)  ?n]  -  wro(n+lKn  (74) 
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when  tp  from  (77)  is  substituted  into  (71)  the  result  can  be  written 
as 


9u/u 


{2<sReL 


— )w  -  <3G7G8+G7) 


+  ( 


^->4 


e,s 


[(- 


3u/u 


16G11G9 


e*  ,  1  .  .  „„  , 2 


'wVu'  *J,+4S«1 

e ,  s  L 


8G9G5 


[(• 


3u/u 


+  G, 


9?  W(62u/  ReT) 

e,s  L 


+4  G13] 


=  0 


(78) 


3u/u 

Equation  (78)  is  solved  for  ( - §.) 


r8u/ue 


3c  'w  <$2ui  sRe, 


by  iteration. 


Je,s^L 

When  a  value  is  taken  for  ( - ,  6  is  known.  In  the  present 

Ar  v’ 

3u/ue  ^ 

investigation  — )w  is  ta^en  as  unity  which  is  consistent  with 

laminar  flow  at  the  stagnation  point  and  a  "power"  profile.  When 
1 

is  known,  a  substitution  into  (77)  gives 


3u/up 
(— -— — — ) 


3?  w  ^ 2  ue , sReL 

Thus  6S  and  *PS  can  be  found  for  a  =  0,  both  for  spin  and  no  spin. 

The  value  of  5S  for  a  ^  0  ,  w  =  0  is  found  by  similar  procedure. 
For  w  =  0,  the  parameters  B  and  C  are  zero  and  A  =  qe.  Near  the 
stagnation  point  the  flow  outside  the  boundary  layer  can  be  expressed 
as  (p.  462  Reference  10) 


u 

e 


clSx 


(79) 


v 

e 


c2se 


(80) 
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where  sx  and  se  are  distances  measured  from  the  stagnation  point 
along  0  =  constant  and  x  =  constant,  respectively.  The  equation  for 
a  streamline  is 


dse 


(81) 


or  with  (79)  and  (80), 


ds  „ 

e 

dsx 


(82) 


where  C3  =  C2/c±.  Equation  (82)  results  in 


s 


0 


c  .s 
4  x 


(83) 


or 


ds 


0 


C.C-S 

4  3  x 


'3-1 


(84) 


At  a  stagnation  point  that  is  off  the  axis  of  revolution  of  an 
elongated  body  of  revolution,  the  radius  of  curvature  of  the  surface 
in  a  plane  9  =  constant  is  greater  than  the  radius  of  curvature  in  a 
plane  x  =  c.  Consequently  c2>ci  and  so  C3>1.  Therefore  all  the 


streamlines  except  the  one  for  04  =  -*>  have 


=  0  at  the  stagnation 

dsx 


point  (see  Figure  4).  Thus  for  all  the  streamlines  except  the  one 
for  04  =  00,  the  angle  a  is  zero  at  the  stagnation  point.  The 
stagnation  point  is  a  nodal  point  of  attachment  (p  76  Reference  10). 
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The  values  for  the  J's  at  the  stagnation  point  are  obtained  by 
putting  (0  =  0  in  (66).  Moreover,  going  to  the  stagnation  point  but 
not  putting  ue  =  0  yet,  it  can  be  shown  that 


3J3 

3x 


=  0 


3J( 

36 


2_  3a 
-  U  G,  ,  7T7T 
e  11  36 


3J, 

_ t 

3x 


-U  2g.  |S 

e  1  3x 


3J. 

TT 


=  2u 


3u 

ep 

e36~^4 


3J-  3  u 

Jx~  =  2uaJx~  G8G7 


3J{ 

W 


2_  3a 
ue  G1236 


!^6  _  _  2  3a 

3x  e  u113x 


3  J_  3u 

_ 9  _  p  ^  /-• 

36  “  e  36  ^11 


(85) 


3  J 


14 


3x 


3ue 

_2u  G, 
e3x  13 


3J 


12 


36 


“ue  (2G13_G2)36 


^  =  ue2(2G9-G8)|f 


3  J 


13 


38 


=  0 


These  relations  together  with 


w. 


=  iIjt 


w. 


are  used  (61)  and  the  resulting  equation  divided  through  by  ue2. 
The  result  is 
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'rocos+MG13  *'rocos*ll5Gl3+ll'l'2su+2*|f«S11 


wx  2  3u  3r 

X  2  ^  +  6  ^vT-  3x  Gi3^r0C0s^+2^Gi333Pc°s4>  +  ^2G11|frocos(})  (86) 


3  cl  o  v  -  u 

+  (2G9-G8>Hroc°s*+;nr  ;rVoCOS+-  *¥  3eSG1i+'-K2GX3-G2)|f 


G„  3 v 

+*  vT  JT  1 

e 


Now  use 


3u/u 

T  ( ~J=-) 

wx  1  3$  ;w 

„  2  ReT  -S  u 

u_  L  e 


(See  Equation  (69)) 


Also,  on  the  line  G  =  tt 


(Appendix  A) 
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and  at  the  stagnation  point 


92u 


(JL?)  = 

v  9  6 

s 


1 

*2 


90' 


9v 

( — — ) 
96  s 


(Appendix  A) 


(89) 


9  2u 


Because  neither  ( 


-) 


nor 


term  (|f)s 


3ve 
(7e_)  s 


is  zero  (Appendix  A)  the 


and  L 
ue 

becomes 


9  v 


=  o 

9x 


9  02  ' s 

is  neither  zero  nor  infinite.  When  (87),  (88),  (89), 
(Appendix  A)  are  used,  the  right  hand  side  of  (86) 


9u/u 
r - - 

1  35 


w 


o  Re_  6  u 
L  e 


0  9u  9ro 

+  6  ltT  3TT  Gi3rooos*  + 
e 


_  .  9u 
2\p  en 

U  90^11 
e 


9a  G2  9ve 

+  (2G13~G2)^g-  +  —  gfg-] 

e 


(90) 


To  prevent  (90)  from  becoming  infinite  as  ue  approaches  zero  at  the 
stagnation  point  it  is  necessary  that  $  =  0  at  the  stagnation  point. 
Therefore  the  value  of  f  at  a  stagnation  point  for  a  f  0,  to  =  0 
is  zero. 

When  ip  =  0  and  (66)  with  w  =  0  and  (85)  are  used  in  (60)  the 
result  is 


roc°s<j)||ue2G7G8+||  6ue' 


G, .=r  t 
14  o  w 


5  (2u 


9u 
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rocos^  G7G8 


2  9r  9u 

Ue  G7G8COS<^— 9x  +  roCOS<^Ue-Tx  G7+Ue  Gl2ft> 


(91) 
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Equation  (91)  is  divided  through  by  u02.  The  term  xw  becomes 

1 


lW' 


U, 


which  is  written  as 


ue5ReL 


3u/ue 

r — - — )  The  relation  u~  =  ui  e-s,, 

'  9?  ' w  e  e's  x 


where  sx  is  measured  from  the  stagnation  point  is  also  used,  as  is 
the  relation  =  GyGg .  Then  (91)  becomes 


-  itlL-r  i  ^  zr 
-o9i^l2  Je  6G 


14 


u 


r  9u/u  r 

— - ( — 3r^)w  -SI—^G  G  +G  > 

s  SReT  3?  w  sx  7  8  7 


e,s  x 


+  g7g8(9s  +  ae)] 


(92) 


In  order  to  prevent  the  right  hand  side  of  (92)  from  becoming  infinite 
as  sx  becomes  zero  it  is  necessary  that 


or 


u  „  6  Re¬ 
ef  o  L 


9u/u 

r  t 

'  9? 


-)w  -  6(2G?G8+G7)  =  0 


(93) 


6s  * 


9u/u 


1 


1 

2 


u  ReT 
e,s  L 


(2G7Gg+G7) 


(94) 


Thus,  for  a  ^  0,  w  =  0  the  value  of  <5  at  the  stagnation  point  is 
given  by  (94).  The  starting  value  of  <$  can  now  be  obtained  from 
(65) ,  (78)  ,  and  (94)  . 

For  ot  =  o  the  starting  value  of  '4>  is  gotten  from  (77)  .  For 
a  o,  w  -■  0  the  starting  value  of  $  is  zero.  For  &  ?  0,  “  ^  o  the 
value  of  ^  at  the  stagnation  point  is  obtained  by  use  of  (39).  Thus 
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for  all  streamlines  except  the  one  for  C4  =  00  (See  (8  3)  and  Figure  4) 
Vg  lies  in  the  plane  x  =  constant  through  the  stagnation  point  because 


a  =  0  there.  Therefore  at  the  surface  Vq 


wro  and  (— g£)<0  (see  (46), 


(47)).  Also  Uq  =  0  for  all  y  at  the  stagnation  point;  therefore 
3uq 

— gy  -  0  for  all  y.  Thus  it  follows  from  (39)  that  'P  =  -  00  for  all 


streamlines  at  the  stagnation  point  except  the  one  for  c 4  =  °°.  For 
04  =  00 ,  a  =  n/2  and  similar  reasoning  results  in  <J>  =  0.  Thus  <j>  is 
double  valued  at  the  stagnation  point.  The  integrations  are  made  for 
C4  =  °°  so  that  ip  =  -  00  is  the  starting  value. 


Friction  Coefficient 

The  component  of  the  friction  coefficient  in  the  direction  of 


qe,  the  velocity  at  the  outer  edge  of  the  boundary  layer,  is 
approximated  by  a  formula  based  on  6  and  on  the  magnitude  of  the 
component  in  the  direction  of  qe  of  the  relative  velocity  of  qe  with 
respect  to  the  surface.  The  formula  is 


Equation  (96)  is  consistent  with  the  expression  (32)  for  the  velocity 
profile  in  the  direction  of  qe.  The  component  tw  along  a  meridian 
is  given  by  (see  Figure  5) 


(cos  a  -  ip  sin 


a) 


(97) 
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and  the  component;  in  a  plane  x  =  constant  by 


(sin  a  +  ip  cos  a) 


(98) 


Integration  of  Equation  for  Displacement  Surface 

The  quantity,  Q,  in  Equation  (26)  is  the  right  hand  side  of 

Equation  (24).  The  integrals  fj*(ue-u)dy  and  /$(ve-v)dy  are 
given  in  Table  II.  The  Equation  (26)  becomes 


d<$* 

dx 
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r  cosbu 
o  Y  e 
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cos  $ 
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3u 
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-0  rr~~  +  - - *-  -r— —  + 
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u  3x 
e 
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o  e 
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3 J,  n  3J 


19.  .  1  36 


(JioTq-  +  ^  ^  +  ITT;  '  TS-(Jio  J,Q) 


r  u  '1830  T 36 
o  e 


30 


r  u  30  18  T  '  “19' 

o  e 


(99) 


The  quantities  6,  ip ,  and  their  derivatives  are  obtained  by  use  of 
(60),  (61),  (62)  and  (63).  The  derivatives  of  the  J’s  are  obtained 


by  calculating  the  J's  at  x+Ax  and  0+A°  and  using  ItZ.  =  —  and 
3J  ~  AJ  -  3r' 


36  A0 


The  velocity  derivatives  and  — —  are  obtained  from  the 


velocity  distribution  and  shape  of  the  body. 
6S 

6*  is  taken  as  —  which  is  valid  for  a  power 


The  initial  value  of 
profile  with  n=l,  the 
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approximation  to  a  laminar  profile  at  the  stagnation  point.  In  (99), 

the  derivative  is  along  a  streamline.  Thus  (27),  (60),  (61), 
dx 

and  (99)  are  integrated  simultaneously  step-by-step.  The  result  is  5, 
\p,  6*,  and  the  streamline  path  0{x).  Different  streamlines  are  ob¬ 
tained  by  beginning  the  integrations  at  points  on  a  circle  of  small 
radius  in  the  tangent  plane  at  the  stagnation  point.  Through  each 
point  on  the  circle  goes  only  one  streamline;  the  entire  streamline 
is  thus  fixed  by  choosing  the  point  on  the  tangent  circle.  The 
relation  for  the  initial  value  of  0  in  the  integration  is 


0 


IT  +  A  0 


IT 


I  . 

-  sm  n 


(100) 


and  the  relation  for  the  initial  value  of  x  is 


x  =  x  +  Ax  =  x  +1  cos<{>  cos  n 
s  s 


(101) 


The  angle  n  is  zero  along  the  meridian  0=u  in  the  direction  of 
positive  x.  The  distance  I  is  the  radius  of  the  startir.  circle. 

Derivation  of  Equations  for  Magnus  Force  and  Moment 

After  the  displacement  surface  thickness  6*  has  been  computed 
it  is  added  to  the  body.  To  calculate  the  force  and  moment  on  the 
resulting  body,  slender-body  theory  (References  5  and  6)  can  be  used. 
In  slender-body  theory  the  force  acting  on  the  portion  of  the  body 
between  the  nose  and  a  section  (C-C)  perpendicular  to  the  £  axis 
(Figure  6)  is  given  by 


(page  50,  Reference  6) 


(102) 
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Here  <J>Q  is  the  perturbation  potential  in  the  zQ  plane  and  is  a 
solution  of  Laplace's  equation. 

Equation  (102)  can  be  placed  in  a  form  more  convenient  for  use 
in  the  present  analysis  by  integrating  by  parts.  Then  (102)  becomes 


~2i  Iz  (A*)  - 
°1  ° 


(103) 


where  A<J>0  is  the  change  in  <i>0  in  going  around  the  section  of  the  body, 
starting  the  ending  at  the  point  7,\.  Because  A4>0  is  zero,  (See 
9.3.15  of  Reference  5),  Equation  (103)  becomes 


F 


z2 


2i 


c 


(104) 


In  order  to  use  (104)  ,  Vg  must  be  known  along  the  contour  c. 
To  find  Vg  along  c  the  velocity  vector  6  is  written  as 


Q  = 


,  i  .  34>.  2 
(x  9^)el 


v  ~—v 

°°o  V 


(105) 


where  4>  is  the  perturbation  potential,  is  a  unit  vector  along  the 
C  axis  which  is  along  and  the  unit  vector  0  is  normal^to  the  body 
cross  section  in  a  plane  t,=c  (Figure  6)  .  The  unit  normal  n  to  the 
surface  of  the  body  is 


=  T 


(106) 
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Then  (108)  becomes 


<m>T5 + $<Thhc 


h  fV 

ho  ds5 


(112) 


According  to  slender-body  theory  the  body  has  a  small  slope  in  a 
plane  B=c,  therefore  T^<<1.  Alsr?  the  disturbance  velocity  ~  is 
small.  Therefore  the  product  *-=■  ^  In  (112)  is  neglected  with  the 
result  that 


Also 


Thhol§^ 


_i  ° 

h_  ds,. 


1113) 


ds^  = 


i  dn^  * 

1+(E-3r»  hodB  (5=c) 


(114) 


T(S,hQ,B)  =  hQ  -  h0U,B)  =  c 
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Then 


3h 

Ts  ‘ 


and 
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3h 

-  ( — 2.) 

^33  ' 
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and 
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1 


Then  (113)  becomes 


3$o 

3v 


3ho 


i  3h„  2 

1+i2(38£) 


(115) 


where  $  in  (113)  has  been  replaced  by  <I>0.  According  to  the  approxi¬ 
mations  of  slender-body  theory,  it  is  permissible  to  replace  the 
three-dimensional  disturbance  potential  $  in  (113)  by  the  two- 
dimensional  disturbance  potential  that  is  a  solution  of  Laplace's 
equation  in  a  plane  £=c. 


When  the  displacement  surface  has  been  added  to  the  body, 

3  4 > 

Equation  (115)  is  used  to  obtain  _ _,  which  equals  Vn  the  non-dimen- 

3v 

sional  velocity  normal  to  a  section  £=c.  The  expressions  for 
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9ho 

( 9  ^  ^  / 


3ho 

( — —)  are  derived  in  Appendix  B.  Equation  (115) 
3  8  £ 


gives  the  non-dimensional  velocity  normal  to  a  body  section  £=c.  In 


the  calculation  of 


3hn  dhn 

, — and  ( _)  the  quantity  £  occurs  (see  B-21 

Ur  fa  '3$  1  ^ 


L3£  8 

and  B-22) .  The  boundary  layer  calculations  provide  6*  and  other 
quantities  as  functions  of  x.  Consequently  £  is  found  from  a  specifi¬ 
ed  value  of  x.  The  value  of  £  for  a  fixed  value  of  x,  however, 
depends  on  6  (See  B-39) .  Consequently  when  x  is  fixed  the  value  of 

34> 

£  varies  slightly  with  6.  The  normal  velocity  therefore  is  ob¬ 
tained  on  a  curve  that  does  not  lie  exactly  in  a  plane  £=c.  By 
putting  R=r0+s*  in  (B-39)  and  noting  that  <s*,  sin  a,  and  r„  are  all 
sma) 1  of  o(a)  ,  it  follows  that  the  products  rQ  sin  a  and  6*  sin  a  are 
of  the  o(a2).  Therefore,  the  departure  from  the  plane 


£c  =  (d+x0)cos  a 
34> 

of  the  curve  on  which  7—^  is  calculated  is  of  o(az),  a  small  quantity 
when  a  is  small.  oV 

If  the  value  of  £  were  fixed  at  the  value  given  by 


£c  =  (d+xQ)cos  a 


the  values  of  x  such  that  £=£c  could  be  found  by  putting 
£=£c=  (dfxQ) cos  a  in  (B-39)  and  solving  for  x  to  obtain 


x  =  xq-R (x, 0 ) cos 0 tana 


(116) 


'X  P  £  JX  T  £ 

The  values  of  -2——  and  tttt-  for  x(0)  would  then  be  found  by  interpola- 


3x 


30 


36’ 


tion  from  the  computed  values  of  * —  and 

o  X 


36* 

30 


at  two  values  of  x,  one 


larger  and  one  smaller  than  the  largest  and  smallest  value  of  x 
given  by  (116).  This  procedure  was  not  carried  out  in  the  example 
calculation  for  the  half-ellipsoid  because  the  additional  complication 
it  would  have  introduced  into  the  calculations  was  judged  not  to  be 


justifiable.  In  the  calculation  of  ( 


3hn  3h 

— -)  ,  ( — -)  etc* 

33  3£ 


(Appendix  B)  , 
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the  exact  expressions  involve  hardly  any  more  computation  than  expres¬ 
sions  that  would  be  obtained  by  neglecting  terms  of  o(az).  Consequen¬ 
tly  the  exact  expressions  were  used. 

The  next  step  in  the  process  of  finding  for  use  in  (104)  is 
to  move  the  origin  of  coordinates  in  a  plane  £=c  back  to  the  axis 
of  revolution  ,  that  is,  to  the  x  axis  (see  Figure  7).  The  equations 
for  doing  this  are 


r  =  +  2hQ(l+x)sin  a  cos  8+  (d+x)  2s.in2a  (117) 


and 


-h  sin  B 

cos  y=  -  (118) 

r 


At  this  stage  in  the  calculation  the  vertical  velocity  along 
the  distorted  circle  with  its  center  on  the  x  axis  and  lying  in  a 
£  =  constant  plane  is  known.  To  find  Vg  from  V'  conformal  transfor¬ 
mation  theory  as  applied  in  Reference  11,  12  and  13  is  used.  The 
method  of  Reference  11  and  12  is  used  to  transform  the  distorted 
circle  in  the  r',  y  plane  into  a  circle  in  the  A,  o  plane  (Figure  7). 
Because 


v'  =  2*0.  -  w 

n  9v  Av 


(119) 


and 


..  O'I'o  /n.  A$o 

n  3 A  ~  AA 


(120) 


and  because  the  potential  at  corresponding  points  in  the  Z'  and 
Z  planes  (Figure  7)  is  unchanged  by  a  conformal  transformation  of 
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the  Z'  plane  into  the  Z  plane  (p  163  of  Reference  14),  it  follows 
from  (119)  and  (120)  that 


But  Av  =  |dZ*  |  and  AA  =  |dz|  ,*  therefore  (121)  is 


(121) 


(122) 


The  ratio  for  Doints  on  the  cross  sections  in  the  Z  and  Z' 

aZ 

planes  is  given  by  (Equation  (37)  of  Reference  12),  namely 


e  o 


1  4-  I  > 

1+  (av} 


(123) 


wnere  for  points  on  the  cross  sections  in  the  Z'  and  Z  planes. 


'  /  iY  Y 

Z  =  r  e  =  r  e 


(124) 


,  io  ftn+io 

Z  =  Ae  =  re  0 
o 


(125) 
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and 


e  =  a  -  Y 


(126) 


The  constant  ft0  is  calculated  by 


(Page  8  of  Refer¬ 
ence  11  and 
Equation  C  of 
Reference  12) 


2tt 


^  / 


ft  (y)dy 


(127) 


The  quantity  e  is  calculated  from 


2  n 


e  (cr  )  «  -■ 


1 

2tt 


(Equation  13 
Reference  12) 


ft  (a)  cot  — 4?-do 


(128) 


Because  6*  is  a  small  fraction  of  rQ,  the  cross  section  in  the  plane 
?  =  c  is  almost  a  circle.  Therefore  the  quantity  o  in  (128)  can  be 
replaced  by  y  (see  Reference  11)  so  that  the  equation  used  to 
calculate  e  is 


e  (y  )  =  -- 


ft(r)  cot  ^'.jy 


(129) 


For  more  distorted  cross  sections  iteration  is  used  (see  Reference  12). 

Once  Vn  is  known  on  the  circle  in  the  A,  o  plane,  the  velocity 
Vs  can  be  calculated  from  an  equation  given  on  page  5  of  Reference  13; 
the  equation  is 
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V°  >  = 


271 

-27  J 


(3  30) 


When  Vs (a)  is  known  from  (130)/  Vs  (r',y)  is  found  from 


v/  ,r'-v)  =  V  (o)lSlb- 1 


(131) 


where  is  the  reciprocal  of  given  by  (123)  .  Then  because 

vs  (t'/Y)  =  Vg(h0f3)  at  corresponding  points,  the  integral  in  (104) 
can  be  evaluated. 

Actually  it  is  more  convenient  to  express  the  integral  in  (104) 
in  terms  of  quantities  in  the  (A, a)  plane.  From  (131) 


va  (r'  '*>  "  Vo) 


(132) 


But  Vs(h0,6)  =  Vs  (r',Y)  and  ds^  =  ds£  because  the  zQ  and  Z'  planes 
are  identical  except  for  a  shift  in  origin  (see  Figure  7).  Therefore 
by  using  (132)  the  term  Vs  ds^  in  (104)  can  be  written  as 


V  dsr  =  V  (a)ds 

S  c,  S 


(133) 


Also,  from  Figures  7  and  8 


z  =  Z  -  (d+x)sin  a 
o 


(134) 
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or  with  r  e 
o 


ft 


=  r/(y) 


and  ds  =  Ado 


(138) 


and 


(139) 


where 


(140) 


and  and  y  are  known  functions  of  a  because  fi=n(y)  and  y =o-c. 

The  force  Fy  is  the  Magnus  force  for  the  section  of  the  body 
between  the  nose  and  the  section  x.  The  force  Fz  lies  in  the  plane 
through  the  free  stream  velocity  vector  and  the  axis  of  revolution 
of  the  body. 
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The  Magnus  moment  about  the  nose  of  the  body  is  given  by  the 
expression,  with  L  taken  as  the  body  length, 


(141) 


where 


|v«,2l2 


is  given  by  (139).  By  integration  by  parts 


Equation  (141)  can  be  placed  in  the  more  convenient  form 


M. 


W 


=  dl 


W 

2  00 


x=0 


(142) 


Once  the  boundary  layer  displacement  surface  surrounding  the 
body  has  been  calculated,  the  Magnus  force  and  moment  can  be 
calculated  by  a  basically  more  exact  method  than  slender-body  theory, 
namely  potential  flow  theory  as  applied  in  Reference  7.  Although 
this  method  was  not  used  in  the  present  investigation  because  of  a 
lack  of  funds,  its  use  is  discussed  in  Appendix  C. 

CALCULATION  METHOD  APPLIED  TO  EQJENESS  RATIO  5  HALF -ELLIPSOID 

To  begin  the  numerical  integration  the  value  of  I  in  TlOO)  and 
(101)  is  taken  as  10-s  for  all  ri  except  ±95°.  For  n=±95°  I  is  taken 
as  10~3  in  order  to  increase  the  spacing  between  the  n=±90°  and 
n=±95°  streamlines  (Figure  9). 

The  increments  in  e  and  6*  are  calculated  from  (27)  and  (99) 
by  use  of  .the  relations 
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In  no  case  is  Ax  allowed  to  exceed  10"4. 

To  calculate  the  derivatives  of  the  J's  (Equation  60,  61)  A0  and 
Ax  are  both  taken  as  10“ 4  and  the  derivatives  are  calculated  by 

l 

Ji,(x+Ax)  Ji,x 
Ax 


Ji , (9+A6 ) ~Ji , 9 
A0 


When  Equations  (27)  ,  (60)  ,  (61)  and  (99)  were  integrated,  it 
was  found  that  rapidly  diverging  oscillations  occurred  at  the  start. 
This  caused  the  computer  to  stop.  Rather  than  spend  time  trying  to 
eliminate  the  divergence,  the  values  of  6,  ip,  and  6*  were  fixed  at 
their  stagnation  point  values  until  x  became  larger  than  -.998.  They 
were  then  allowed  to  vary  and  no  oscillations  occurred.  The  stag¬ 
nation  point  is  at  -.999910  so  that  6,  \p,  and  6*  were  fixed  for  the 
first  .001910  in  x.  For  the  largest  value  of  n ,  namely  ±95°,  a 
calculation  shows  that  the  distance  along  the  streamline  is  about 
.0057  compared  with  .00191  along  x.  Thus  for  the  streamlines  calcu¬ 
lated  (Figure  9),  6,  ip,  and  S*  are  fixed  for  distances  along  stream¬ 
lines  not  greater  than  .0057  from  the  stagnation  point.  Consequently 
the  results  for  6*,  and  so  for  the  Magnus  force  on  the  portion  of  the 
body  near  the  nose,  may  be  inaccurate.  Effects  of  changes  in  6,  \p, 
and  6*  near  the  stagnation  point  decrease  rapidly  with  increase  in  x. 
Thus,  a  change  in  6  at  the  stagnation  point  by  a  factor  of  5  for 
n=95°  makes  only  a  change  of  3  in  the  fourth  significant  figure  in 
6*  for  x  =  -.600. 


3J. 

3x 
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To  evaluate  the  integral  in  (129)  and  (130)  the  method  of 
Reference  11  is  followed.  Equation  (129)  is  integrated  by  first 
noting  that  the  integrand  becomes  infinite  at  y'  =  y.  Consequently 
(129)  is  written  as 


e  (y') 


JL 

2tt 


/ 

Q  (y)  cot  -Y-^-dy 


+ 


+ 


1 


i 


i 

i 
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r  /  y«no  ,  y2 

+  /  <  /  ^(Y)cot^-dY  +|  ft(y) 

0o^°  /  i t  *V+n0 


cot  — dy 


Then  near  y=y* ,  ft(y)  is  written  as 


fl(Y)  =  A(y')  +  (Y-Y7)  f^y')  +  •  • 


and  substituted  into 


o  y0 

/  /■  2  / 

fl  (Y )  cot^'  dy+  |  n  (Y )  cot dy 


y/f  n. 


dO  /  / 

thereby  giving  4 A  g^-(Y  )  /  where  y  -y  j  =  A=y2“Y 


Then  equation  (143)  becomes 


e  (y  )=  - 


1 

2tt 


/ 

y  -A 
J"  fl(Y) 


ft  (y)  cot^"-|— dy  + 


+  4Oy 


'■] 


2n 

/ 

ft  (y) cot^-^-dy 

/ 

y  +A 


(143) 


(144) 


49 


NOLTR  72-80 


50 


NOLTR  72-80 


where  N=—  .  For  n=0,  instead  of  (146)  the  result  is 
A 


e  <°>  ■  ‘2? 


N-l 


4A  |f< 


k=2 


Q,  In 
k 


o-  kA 
Sm  2 


Sin 


(k-l)A 


(147) 


Also  e(N)=e(o).  In  the  calculations,  N  was  taken  as  180.  A  test 
with  N=360  showed  a  noticeable  but  negligible  change  in  e  (y) . 

To  evaluate  the  integrals  in  (138)  and  (139)  the  trapezoidal 
rule  was  used  with  the  interval  2tt  divided  into  180  equal  parts. 

The  program  for  the  numerical  calculations  for  use  on  the 
CDC  6400  is  given  in  Table  IV. 

RESULTS  FOR  FINENESS  RATIO  5  HALF -ELLIPSOID 

As  an  example,  the  Magnus  force  and  moment  on  a  fineness  ratio 
5  half-ellipsoid  spinning  at  p=.25  and  at  a=4°  is  calculated.  The 
boundary  layer  is  turbulent  and  is  characterized  by  n=l/5.  The 
Reynolds  number  based  on  the  length  of  the  half-ellipsoid  is 
4. 8x10 6 . 

The  inviscid  flow  velocity  distribution  over  the  half-ellipsoid 
is  assumed  to  be  the  same  as  over  the  forward  half  of  a  complete 
ellipsoid.  Because  of  the  wake  this  assumption  is  not  exact, 
especially  near  the  base,  but  is  used  in  order  to  obtain  an  analytic 
expression  for  the  velocity  distribution  in  the  example  calculation 
which  is  of  an  exploratory  nature.  The  velocity  distribution  is 
obtained  from  Chapter  V  of  Reference  15.  The  paramete  1  (Figure  8) 
equals  unity  here.  The  result  is 


u  = 
.  e 


1— x* 


2  2 
l-xz(l-tZ) 


1 

2 


cos  a 


[l-x2(l-t2)]2 


tK2  cos  6  sin  a  (148) 
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and 


v  =  -K»  sin  0  sin  a  (149) 

e  2 


where 


and 


The  half-ellipsoid  for  which  the  calculations  are  made  has  t=.l  so 
that  Ki=l. 02070  and  K2=l. 96023.  The  shape  of  the  half-ellipsoid  is 
given  by 


Calculated  streamlines  on  the  half-ellipsoid  are  shown  in 
Figure  9  for  positive  values  of  n .  Each  streamline  shown  in  Figure  9 
has  a  corresponding  streamline  which  is  its  reflection  in  the  line 
0=tt  and  which  has  a  negative  value  of  n. 

The  boundary  layer  displacement  thickness  6*  along  the  stream¬ 
lines  shown  in  Figure  9  (n>o)  are  presented  for  n>o  in  Figure  10a 
and  for  n<o  in  Figure  10b. 
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In  Figure  11a  and  lib  is  given  the  distribution  of  the  full 
boundary  layer  thickness  6  along  various  streamlines. 

In  Figure  12a  and  12b  is  shown  the  variation  of  5*  with  0  at 
x=-.5  and  at  x=0.  The  variation  is  slightly  unsymmetric  about 
0=rr ;  it  is  this  asymmetry  that  produces  the  Magnus  force. 

In  Figure  13,  (a)-(i)  is  shown  the  variation  of  tp  with  x  for 
all  the  streamlines  calculated.  The  calculations  were  begun  with 
tp=- 1  at  the  stagnation  point  instead  of  tp=-°°.  A  few  trials  showed 
that  in  the  present  example,  a=4°  and  p=.25,  when  the  initial 
value  of  ip  was  less  than  -5  or  so,  the  integration  of  the  differential 
equations  for  6*,  6,  and  tp  could  not  be  made  because  of  divergent 
oscillations.  The  behavior  of  tp  near  the  stagnation  point  shown  in 
Figure  13  indicates  that  the  effect  of  the  choice  of  the  initial 
value  of  ip  soon  disappears.  Moreover,  two  calculations  for  a=6°, 
p=l,  n=-10,  one  with  tps=-5  and  the  other  with  t|>s=-10  showed  that 
when  x  had  increased  to  -.985555  from  the  stagnation  point  value 
-.999796,  the  values  of  tp  for  the  two  calculations  differed  by  only 
1  in  the  third  significant  figure.  The  difference  then  decreased 
as  x  increased. 


In  Figure  14  is  shown  the  normal  force  coefficient  based  on 
the  local  cross  section  area  for  the  portion  of  the  body  between 


station  x  and  the  nose.  The  coefficient  Cn  is  obtained  from  - 

(Equation  (137))  by  the  relation  Pa>Vco2-? 

_  l 


(152) 


In  Figure  15  is  shown  the  corresponding  Magnus  force  coefficient  Cy 

Fy 

obtained  from  — — -  (Equation  (138))  by 

PooVoo  ^r2 

—  Xj 


(153) 
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Although  the  sign  of  Cy  computed  in  the  present  method  is  positive 
along  the  +Y  axis  in  Figure  6  for  w  as  shown  in  Figure  1,  it  is  shown 
negative  in  Figure  15  in  order  to  adhere  to  the  sign  convention 
given  in  Reference  16. 

The  curves  for  and  Cy  are  shown  dashed  for  x<-.6  because  from 
calculations  for  n-±95°  and  a  calculation  for  n=+95°  with  6S  equal 
to  5  times  the  correct  value  of  6S,  it  was  found  that  the  change  in 
6*  for  x<-.6G0  caused  by  using  56 s  instead  of  6S  in  the  calculation 
for  n=+95°  was  greater  than  10  percent  of  the  difference  between  6* 
for  n=95°  and  6*  for  n=~95°.  Thus  the  Magnus  force  for  x<-„600  may 
be  affected  by  more  than  10  percent  by  inaccuracies  in  the  initial 
values  of  6  and  ip  and  by  the  method  of  numerical  integration;  see 
section  on  "Numerical  Calculations". 

The  results  in  Figures  (14)  and  (15)  indicate  that  because  of 
the  boundary  layer,  C$j  decreases  and  Cy  increases  as  the  fineness 
ratio  increases.  The  decrease  in  Cn  is  almost  linear  with  x  for 
larger  x  but  Cy  increases  more  rapidly  than  in  a  linear  manner  all 
the  way  to  x=0,  the  body  base. 

The  value  of  the  moment  coefficient  of  the  entire  body  defined 
by  (Reference  16 ) , 


C 


n 


M 


f-V^L3 
2  00 


(154) 


is  found  to  be  1.058xl0“3.  (See  Equation  (142)) 


COMPARISON  WITH  EXPERIMENT 

There  do  not  seem  to  be  any  experimental  data  for  a  spinning 
half-ellipsoid  of  fineness  ratio  5.  Consequently,  strictly  speaking 
there  can  be  no  comparison  of  the  calculated  results  with  experiment. 
There  are,  however,  some  data  for  pointed  shapes  of  near  5  fineness 
ratio  in  subsonic  flow.  One  set  of  data  is  that  of  Reference  16. 

Tests  were  made  of  a  model  with  a  2  cal.  secant-ogive  nose  and  a 
3  cal.  cylinder  after-body.  For  a=4°,  M=.2,  and  ReL=4.8xl06, 
the  data  in  Figure  6  of  Reference  16  indicate  that  Cy  was  about  -.0075. 
This  value  is  about  24  times  as  large  as  the  calculated  value 
-3.11xl0-4.  The  experimental  value  of  Cn  was  about  .148  (Figure  11 
of  Reference  16);  the  calculated  value  is  .1385.  The  calculated  value 
is  93  percent  of  the  experimental  value.  Slender-body  theory, 
without  any  boundary  layer,  p.  67  of  Reference  6,  gives  the  value 
2a  or  .1396.  This  is  almost  equal  to  the  calculated  value  with 
boundary  layer. 
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In  Reference  17  the  results  of  tests  of  a  4,9  caliber  smooth 
body  are  presented.  The  nose  is  an  almost  flat  surface  of  about 
.198  calibers;  this  is  followed  by  a  16-degree  half-angle  cone  for 
about  .775  calibers;  this  in  turn  is  followed  by  a  tangent  ogive 
for  the  next  2.5  calibers.  Tne  rear  portion  is  a  cylinder.  The 
Reynolds  number  based  on  body  length  was  about  3.4xl06;  at  this 
Reynolds  number  the  shape  of  the  nose  portion  almost  certainly 
ensures  turbulent  flow  rearward  from  just  behind  the  nose  cap.  For 
this  body  at  a  Mach  number  of  .8,  the  lowest  tested,  the  value  of 
Cy  was  found  to  be  about  -.0125.  The  change  in  body  shape  and  Mach 
number  from  the  tests  of  Reference  16  to  those  of  Reference  17 
resulted  in  a  change  in  Cy  from  -.0075  to  -.0125.  The  value  computed 
for  the  half-ellipsoid  is  -3.11xl0“4.  Consequently  it  appears  that 
the  calculated  value  of  Cy  for  the  half-ellipsoid  is  probably  much 
too  small.  The  half-ellipsoid  may  have  a  more  favorable  pressure 
distribution  than  the  bodies  tested  and  thus  a  thinner  boundary 
layer.  If  so,  it  would  have  a  smaller  Magnus  force.  It  is  not, 
however,  likely  that  this  is  the  explanation  of  the  large  difference 
between  the  measured  and  calculated  values  of  Cy. 

The  calculated  value  of  Cn,  the  Magnus  yawing  moment  coefficient, 
Equation  (154),  is  1.058xl0-3.  The  value  of  Cn  calculated  by  use  of 
Figure  7  of  Reference  16  is  about  -1.7xl0-3.  The  experimental  values 
become  zero  at  a=6°  and  positive  for  larger  values  of  a.  No  values 
are  shown  for  a<4°.  The  value  of  Cn  tor  the  body  of  Reference  17 
(Figure  42)  is  about  1.61xl0-2. 

DISCUSSION  OF  METHOD 

The  present  method  computes  a  displacement  surface  by  the  use 
of  boundary-layer  theory.  Consequently  all  the  boundary  layer 
approximations  are  inherent  in  the  method.  In  order  for  boundary- 
layer  theory  to  be  applicable  the  boundary- layer  thickness  at  a 
point  on  the  surface  should  be  much  smaller  than  the  distance  from 
the  point  to  the  stagnation  point.  In  the  calculation  for  the  half¬ 
ellipsoid  the  largest  6  occurs  at  x=0,  0=9.  From  Figure  11  it 
appears  that  6  is  about  25  percent  of  rQ  at  x=0,  0=0.  Although  6  is 
a  large  fraction  of  rQ,  it  is  less  than  2.5  percent  of  the  distance 
to  the  stagnation  point  (see  Figure  9).  Therefore  the  boundary  layer 
is  thin  as  required  by  boundary-layer  theory. 

When  the  boundary  layer  approximations  are  made, it  can  be  shown 

K-,  t 


3J? 

«y 


•cos 


<p 


ap 

When  v  is  of  the  same  order  as  u  (u=o(l))  ,  ■gy  is  of  order  unity 
instead  of  o(S)  as  is  usual  without  spin  (6<<1) .  Therefore  the 
static  pressure  difference  across  the  boundary  layer,  which  is 

proportional  to  can  be  of  o(<5)  for  a  body  with  spin  instead  of 
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o(62)  as  is  usual  without  spin.  Although  the  difference  in  static 
pressure  across  the  boundary  layer  is  of  o(6)  instead  of  o(S2),  it 
is  neglected  in  the  present  analysis  because  terras  of  o  (6)  have  been 
neglected  in  deriving  the  boundary  layer  equations  from  the  complete 
equations  of  motion  and  because  its  inclusion  would  have  made  the 
method  more  complicated.  Kelley  and  Thacker  in  their  analysis  for  a 
cylinder  with  a  laminar  boundary  layer  find  that  the  effect  of  the 
pressure  gradient  across  the  boundary  layer  changes  the  Magnus  force 
by  about  30  percent.  Although  this  is  appreciable,  for  the  half¬ 
ellipsoid  an  effect  of  this  amount  would  not  change  the  impression 
that  the  calculated  value  of  the  Magnus  force  coefficient  is  much 
too  small. 

In  the  present  metnod,  the  velocity  profile  shape  parameter  n  is 
fixed  before  a  calculation  is  made  and  does  not  vary  along  a  stream¬ 
line.  If  n  were  allowed  to  change  by  introducing  an  additional 
equation  for  the  change  of  n,  the  displacement  thicknes  5*  would 
respond  more  accurately  to  pressure  gradients  along  the  streamlines. 
The  introduction  of  an  equation  for  the  change  in  n  would,  however, 
further  complicate  an  already  complicated  method  and  introduce 
additional  uncertainties  associated  with  information  needed  for  the 
n  equation. 

The  assumption  that  the  u  and  v  velocity  profiles  have  the  same 
6  has  been  shown  by  Cooke  (Reference  18)  to  be  a  source  of  inac¬ 
curacy.  The  use  of  two  6's  would  introduce  another  dependent 
variable  and  a  partial  differential  equation  for  it.  The  feeling  is 
that  instead  of  introducing  an  equation  for  n  and  an  equation  for 
another  6  it  might  be  better  to  investigate  the  use  of  a  finite 
difference  method  instead  of  the  integral  approach.  Of  course  this 
introduces  other  problems,  one  of  which  is  the  effect  on  the 
turbulence  of  a  spinning  wall.  Although  this  seems  not  to  appear 
in  the  present  analysis,  it  is  hidden  in  the  friction  formula  and 
velocity  profile  shape,  both  of  which  involve  assumptions. 

Tiie  use  of  slender-body  theory,  a  theory  for  very  thin  bodies 
at  very  small  angle  of  attack,  to  calculate  the  Magnus  force  further 
decreases  the  accuracy  of  the  results.  At  first  glance,  it  appears 
that  more  accurate  results  can  be  obtained  by  using  the  method  of 
Reference  "/  and  8.  Although  slender-body  theory  requires  no  infor¬ 
mation  for  the  wake,  the  method  of  Reference  7  and  8  cannot  be  used 
unless  the  distribution  of  the  displacement  thickness  behind  the 
blunt  base  is  known.  The  method  of  Reference  7  and  8  cannot  be  used 
to  calculate  the  flow  near  the  base  of  a  body  with  a  blunt  base  and 
zero  wake  displacement  thickness;  infinite  velocities  are  obtained 
at  the  corner.  Because  the  behavior  of  the  wake  is  unknown  and 
because  the  Magnus  force  is  small,  the  inaccuracy  caused  by  an 
inexact  representation  of  the  wake  can  be  important.  For  a  body 
with  spin  and  angle  of  attack  the  wake  is  not  axisymmetric.  It 
seems  that  in  order  to  get  an  accurate  solution  it  may  be  necessary 
to  solve  the  entire  flow  problem,  boundary  layer  plus  wake. 
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Another  source  of  error  is  the  lack  of  a  reliable  method  for 
calculating  where  on  the  body  the  initial  laminar  boundary  layer 
changes  to  a  turbulent  one.  Consequently  unless  the  transition 
point  is  fixed  or  known  in  some  way /  the  initial  conditions  for  the 
turbulent  boundary  layer  and  the  location  of  its  origin  are  also 
unknown.  Therefore  even  if  the  entire  flow  field  could  be  solved 
for  exactly,  the  results  might  still  be  unreliable. 

For  tests  of  analytical  methods  for  computing  Magnus  forces  it 
would  be  desirable  to  test  a  simple  shape  such  as  a  half-ellipsoid 
at  a  low  Mach  number.  The  advantage  of  a  shape  like  a  half-ellipsoid 
is  that  the  potential  flow  velocity  distribution  over  at  least  the 
forward  portion  is  given  by  an  analytic  expression.  This  is  a 
help  in  obtaining  velocity  derivatives  and  in  finding  the  flow  near 
the  stagnation  point  at  angle  of  attack. 
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FIG.  I  FLOW  ABOUT  SPINNING  AXISYMMETRIC  BODY  AT  SMALL  ANGLE  OF  ATTACK 
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FIG.  4  STREAMLINES  AT  STAGNATION  POINT 
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Zo(h,0)  Z'(r',y)  Z(X,«r) 


x=C  PLANE  £  =  C  PLANE  £  =  C  PLANE  CIRCLE  PLANE 


FIG.  7  CROSS  SECTIONS  USED  TO  FIND  V's  FOR  USE  IN  SLENDER-BODY  THEORY. 


AZIMUTH  ANGLE  B 


FIG.  9  STREAMLINES  ON  FINENESS  RATIO  5  HALF-ELLIPSOID  AT  4  DEGREES  ANGLE  OF 
ATTACK  FOR  VARIOUS  STREAMLINE  STARTING  ANGLES  ri 


si 

S! 


FIG.  11a 


BOUNDARY  LAYER  THICKNESS  8  ON  HALF-ELLIPSOID  ON  VARIOUS 
STPEAMLINES  SPECIFIED  BY  STARTING  ANGLE  it  (r>>  0). 


FIG.  lib  BOUNDARY  LAYER  THICKNESS  8  ON  HALF-ELI  IPSOID  ON  VARIOUS  STREAMLINES 
SPECIFIED  BY  STARTING  ANGLE  tj  (  t?  <  0). 


TANGENT  OF  ANGLE  BETWEEN  POTENTIAL  FLOW  STREAMLINE  AND  SURFACE  SHEAR  STRESS 
R  STREAMLINES  WITH  STARTING  ANGLES  77  =  i  0.01°. 


FIG.  13f  ty,  TANGENT  OF  ANGLE  BETWEEN  POTENTIAL  FLOW  STREAMLINE  AND  SURFACE  SHEAR  STRESS 
FOR  STREAMLINES  WITH  STARTING  ANGLES  17  =t40°. 


TANGENT  OH  ANGLE  BETWEEN  POTENTIAL  FLOW  STREAMLINE  AND  SURFACE  SHEAR  STRESS 
?  STREAMLINES  WITH  STARTING  ANGLES  ■»?  =  ±  90°. 


FIG.  14  NORMAL  FORCE  COEFFICIENT  CN  BASED  ON  LOCAL  CROSS  SECTION  AREA  FOR 
PORTION  OF  BODY  BETWEEN  STATION  X  AND  THE  BODY  NOSE. 


NOLTR  72-80 


NOLTR  72-80 


APPENDIX  A 


To  obtain  an  expression  for  —  at  the  stagnation  point  for  a  /  0, 
the  expression  3x 


tan  a  = 


(A-l) 


is  used.  By  symmetry,  the  stagnation  point  lies  on  the  curve  9  =  it 
and  ve  =  0  on  6  =  it  for  all  x.  Therefore,  a  =  0  for  9  =  ir  for  all  x 
therefore,  =  0  on  9  =  n  and  therefore  at  the  stagnation  point. 

9v 

Also  e  =  0  for  all  x  so  that  the  term 


-  v  v 

-  - G_r  cos  <}> 

u„  9x  7  o 

e 


(A- 2) 


in  (86)  is  zero  for  9  =  it. 

To  obtain  an  expression  for  —  at  the  stagnation  point,  use  (A-l) 
to  get  3e 


e  90 


2  2 
u  +  v 
e  e 


For  x  = 


x  and  0  near  v, 
s 


3u  32u 

ue  *  <TTT>  „  i0  + 

s  .90 


(A- 3) 
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3v  32v  Aft2 

/  ev  AA  ,  ,  ex  AG 

30  s  3fl2  2 

s 


3u  3u  3  u 

sir  ■  <*r>  +  (— i>  49  +  •  • 

s  30  _ 


3v  3v  3  v 

W~  =  <i~)  +  (“T72)  A0  +  *  * 
s  30 


Because  ue  is  symmetric  with  respect  to  0  =  n  it  follows  that 

3ue  ,  3ue  32ue 

(r~)  =  0  for  0  =  7T ;  therefore,  ( — -)  =  0.  Also  ( - »)  ^  0  for 

30  30  s  90^  s 

a  /  0.  Moreover,  ve  is  antisymmetric  with  respect  to  0  =  it. 

32ve  32ve  3ve 

Therefore,  ( - j)  =  0  for  0  =  tt;  thus  ( - ~)  =  0.  Also  ( - )  f  0 

1  A  /  A  A 


for  a  ji  0.  Then  (A-3)  becomes 


30  s 


3  u  ..,2  3v  3v  3  u 

“<— !>  40 

30  s  s  30 

_ s _ s 

32u  2  4  3v  2 

(~ ~2>  T  +  30-  40 

30  _  s 
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APPENDIX  B 


From  Figure  8  the  vector  S  to  a  point  on  the  body  is  written,  with 
all  lengths  non-dimensionalized  by  the  body  length  L,  as 


=  i  (d+x)  +  Rr 


(B-l) 


also 


exS  +  hQh 


(B-2) 


or,  with  (B-l)  and  (B-2) 


i  (d+x)  +  Rr  =  e^  +  hQh 


(E-3) 


Then 


h  =  i*h(d+x)  +  Rh*r 
o 


(B-4) 


The  unit  vector  h  is  given  by 


h  =  e2  cos  6-  e^  sin  8 


(Figure  8)  (B-5) 


and  the  unit  vector  r  is  given  by 


r  =  j  cos  0  -  e,  sin  0 


(Figure  8)  (B-6) 


B-l 


Then  from  (B-4) ,  (B— 5 ) ,  and  (B-6) 


R  =  R(  £,h,B) 

0  =  d(Z,  h,B) 

Then 

dx  =  x^d$  +  x^dh  +  x^dg 
dR  =  R^d£  +  Rhdh  +  R^dg 

de  =  e^d?  +  ehdh  +  e^dg 

Also 

dH  =  HRdR  +  Hxdx  +  HQd0  *  H^d£  +  Hhdh  +  Hgdg 


(B-10) 

(B-ll) 

(B-12) 

(B-13) 

(B-14) 

(B-15) 


B-2 


w*!?rmq 


H5  ■  V@  +  Vb  +  H96B 


F ix  C5/hQ/3)  =  c  it  follows  that 


(B-18) 


3h 

C_}  = 


Hh 


(B-19) 


and 


8ho  HR 

(33  >  Hh 


(B-20) 


When  (B-16) ,  (B-17) ,  (B-18)  are  used  in  (B-19)  and  (B-20)  the  result 
is 


The  quantities  on  the  R.H.S  of  (B-21)  and  (B—22)  are  obtained 
as  follows.  Equation  (B-8)  is  written  as 


H  ( X / R,  0)  =  R-R  (x,  6 )  =  C 

(B-23) 

Then 

hr  ■  1 

(B-24) 

h  -  (iSi 

Hx-  W# 

(B-25) 

He  “  Ox 

(B-26) 

From  R  =  r  +  6* 
o 

it  follows  that 

d> 

■  ‘Tx1,  +  l2S5>, 

(B-27) 

Also 

d5* 

dx 

.36*.  .36*.  d6 

"  V  3xfl +  1  30' „  dx 

U  X 

(B-28) 

Therefore 

(«i) 

d6*  .36*.  d8 

~  dx  "  1  38;  Hx 

A 

(B-29) 

d6* 

The  quantity  is  computed  during  the  integration  of  the  equation 

35* 

for  the  displacement  surface.  The  quantity  (-jpg—)  is  found  from  the 

x 
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computed  values  of  6*  at  various  6  for  a  fixed  x  by  use  of  the 
spline-fit  method  of  Reference  20.  The  quantity  §|  is  the  stream 
line  direction.  Equation  (B-29)  is  used  with  (B-27)  to  get  (£?-) 

OX  „ 


For  H  e  =  (|~)  there  is  used 


(££)  =  / _ o,  ,  ,36*. 

lae;  +  Hra) 


3r 

where  (j~) 


-  0  because  r  is  independent  of 


Prom  (B-3) ,  (B-S) ,  and  (B-6)  it  follows  that 

R  =  5  sin  a  cos  0  +  hQ (cos  a  cos  6  cos  0  +  sin  0  sin  0)  (B-3 


Then 


R^  =  sin  a  cos  0 


(B-3 


Rh  cos  a  cos  ^  cos  9  +  sin  8  sin  0 


(B-3. 


R8  hQ(-cos  a  sin  8  cos  6  +  cos  8  sin  0) 


Also  from  (B-3),  (B-5)  and  (B-6)  it  follows  that 


(B-3 


x  -  £  cos  a  -  ho  sin  a  cos  8  -  d 


(B-3! 


rhen 


K-  =  cos  a 


(B-3( 


x ,  =  -sin  a  cos  8 


/  rj  n 


x0  =  h  sin  a  sin  0 
p 


(3-38) 


To  find  e^,  0^,  and  0^  multiply  (3-3)  by  e^».  The  result  is 

£  =  (d+x) cos  a  +  R  sin  a  cos  0  (B-39) 

Now  multiply  (B-3)  by  e^*  •  The  result  is 

R  sin  0  =  hQ  sin  0  (B-40) 

By  solving  (B-39)  for  cos  0  and  (B-40)  for  sin  0  there  is  obtained 


h  sin0 

tan  0  =  -= - - - —~r - 5 - 

5  sm  a  +  hQ  cos  0  cos  a 

Equation  (B-3)  is  multiplied  by  e2*  with  the  result 

- (d+x) sin  a  +  R  cos  a  cos  0  =  h  cos  0 

o 

From  (B-40)  and  (B-42) 


(B-41) 


(B-42) 


R  sin  0 


tan  0  = 


-  (d+x)  sin  a  +  R  cos  a  cos  0 


(B-43) 


From  (B-41)  there  is  found 


65  = 


-ti  sin  0  sin  a 
o 


2  2 

U  sin  a  +  hQ  cos  0  cos  a)  +  (hQ  sin  0) 


(B-44) 


B-6 
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£  sin  a  sin  8 

1  =  - 3 - 5 —  (B“45) 

(C  sin  a  +  hQ  cos  8  cos  cr  +  (hQ  sin  8) 


h  (£  sin  a  cos  8  +  h  cos  a) 

0  =  - 2 -  & - —  (b**46) 

p  (£  .sin  a  +  hQ  cos  8  cos  a)  +  (hQ  sin  8) 


k 

£ 
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APPENDIX  C 


Unlike  slender-body  theory  the  method  of  Reference  7  and  8  does 
not  require  the  body  to  be  slender  or  the  angle  of  attack  to  be 
small.  In  contrast  to  slender-body  theory,  however,  the  contour  on 
which  the  pressure  is  computed  cannot  have  a  blunt  base  because  the 
pressure  on  the  body  near  the  base  is  affected  by  the; shape  of the  wake 
behind  it.  A  blunt  base  results  in  infinite  velocities  at  the  base 
corner  and  too  high  velocities  in  the  entire  region  near  the  base. 
Consequently,  the  body  must  be  extended  beyond  the  base  in  a  more 
or  less  realistic  manner.  As  an  approximation  to  the  real  flow  the 
body  can  be  extended  beyond  the  base  by  use  of  the  equation 


(C-l) 


where  =  R^G)  and  the  distance  x  in  (C-l)  is  zero  at  the  body  base 
and  positive  to  the  rear.  The  non-dimensional  radius  R*w  is  the 
wake  displacement  radius  very  far  behind  the  body.  It  depends  on 
the  drag  coefficient  of  the  body;  the  relation  is 


(Appendix  D) 


(C-2) 


The  arbitrariness  in  the  choice  of  N  and  in  equation  (C-l) 

itself  make  it  inadvisable  to  introduce  the  additional  complication 
*  dR*w 

of  R  ^  0  and  (— ~^)  \  0.  Consequently ,  instead  of  (C-l),  there 
can  be  used 


cos 


.2  7 r 


(C-3) 


Calculations  for  tne  body  of  reference  17  by  a  program  of  the 
method  of  Reference  8,  prepared  by  Hess  and  Clissold  of  Douglas 
Aircraft  showed  that  N=20  gave  the  most  satisfactory  velocity  dis¬ 
tribution  over  the  rear  portion  of  the  body. 


When  the  Mach  number  is  not  zero,  the  pressure  coefficient  C 


on  the  body  is  calculated  by  using  Gotherts  Rule  (page  397  of  Ref¬ 
erence  19)  together  with  the  method  of  Reference  8  as  programmed  for 
an  electronic  computer  in  Reference  7.  The  velocity  distribution 
is  calculated  first  for  a  body  of  thickness  ratio  t"||  1-M^^1  at  an 
angle  of  attack  of  a  )J  1-Mm^  The  velocities  are  then  multiplied  by 

— =•  .  This  velocity  distribution  is  used  to  compute  the  boundary 

1-*L“ 

00 

layer  displacement  surface  <5*  which  is  added  to  rQ  to  get  R.  The 


and  the  program  of 


resultant  radius  R  is  then  multiplied  by 
Reference  7  used  to  compute  at  <x| 1-Mto^* ; 

multiplied  by  — ^—5-  to  get  the  C  distribution  used  to  compute  the 
1-M  z  p 


this  C  is  then 
P 


Magnus  force  and  moment. 


The  force  in  the  Y  direction  is  given  by 


C-2 


( 


~  ~  J  P  e3<^A 

where  ft  is  the  outward  unit  normal  to  the  body  surface  plus  the 
displacement  surface  and  e^  is  a  unit  vector  along  the  Y  axis 
Figure  6.  The  unit  vector  n  is  given  by 


n  = 


VH 


Wf 


where  the  equation  of  the  body  is 


H  (x,R, 0)  =  R-R (x,  0 )  =  c 


AAA 


Then  in  the  i,r,0  system  (Figure  8) 


or 


Then 


therefore 


K, 


VH  =  iHx  +  HEr  +  ^  e 


R, 


VH  -iKx  +  r  -  e 


R 


n  = 


-iRx  +  r  -  -t  6 


PI 


2 "  1  +  ("l)2 


R. 


A  A 

n*e. 


-sin  0  +  -=r~  cos  0 

K 


R.2  ♦  1  ♦ 


When  the  static  pressure  P  is  written  as 


where 


dA  =  R 


1  +  (H~cos  *)2 

~  cos 


d0dx 


(C-14) 


has  been  used.  The  non-dimensional  area  dA  is  on  the  body  plus  the 
displacement  surface. 


3  9  3R 

By  using  coscW—  =  -5 — ,  R  =  * — s  =  s  and  cos  <p-+o  near  the  nose  of 

dX  d  S  dSX  X 

X  X 

a  blunt  body  it  can  be  shown  that  at  the  stagnation  point  the 

9  6  *  96* 

integrand  becomes  equal  to  -C  (—  ■)  cos  0 .  The  value  of  (jg— ) 

*s  s  s 

is  zero  because  the  6*  surface  is  assumed  to  be  continuous;  therefore 
the  integrand  in  (C-13)  is  zero  at  the  nose  of  a  blunt  body. 


The  Magnus  moment  about  the  nose  (x=-d)  is  given  by 


sin  0  -  jT^cos  6  ^Jl  +  (|^oos  <{>) 


Poo—  2*"3 

T~Voo  L 
2  00 


Rx2  + 1  +  <4>2 


cos  p 


(C-15) 


-(d+x)d0dx 


where  My  is  positive  as  shown  in  Figure  1. 


APPENDIX  D 


Far  behind  the  body  the  flow  is  parallel  to  V^.  The  wake  dis¬ 
placement  radius  R^*  far  downstream  is  defined  by 


(D-l 


In  (D-l)  y  is  measured  from  the  wake  center  line  and  u  is  parallel 
to 

00 


Equation  (D-l)  can  be  written  as 


R  * 
w 


J  PcoVJdy  -  J  Pjjdy  =  f 

O  o  J 


Pa,v«>ydY  =  /  puydy 


or 


_  R  *2  f 

V  ~ —  =  /  ( p  V  - 

’  00  2  j  4PC0  00 


pu)  ydy 


?hen 


R  *  2 
(~)  = 


00 

r 


(l-pu) 


(D-2) 


2 


Then  (D-6)  becomes 


or 


R  *2 
w 


2  [l+(K-l)Ma2] 


or  with  (D-5) 


R  * 
w 


r„ 


—  1  - n 

[1  +  (H)Mm2] 


MllO 


f  f  i  m  »rW"  7*1 


-  X-* 

5  rM  rM 


(D-ll) 


(D-12) 
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TABLE  I 

Expressions  for  the  Gi(n) 

n-1 _ 

G1  =  Tn+1)  (n+2)  (2n+l) 

1 

G2  =  (n+1)  (n+2) 

n2  +  n-1 
G3  =  (n+1) (n+2) 

n+1 
G^  =  n+2 

3n  +  4 _ 

g5  =  2  (n+1)  (n+2 )  (2n+3) 

-2 (n+1) 2 
Gg  =  Tn+2)  (2n+3) 

n 

G^  =  n+1 
1 

Gg  =  2n+l 
1 

Gg  =  2  (n+1 ) 

-n  (2n+3) 

G2q=  2 (n+1) (n+2) 

-1 

GU=  Tn'+l)  (2n+l)  (2n+3)' 

n _ 

G22=  (n+1) (2n+l) 

1 

G13=  2  (n+1)  (2n+l) 

3n _ 

g14=  ^  (n+l)  (n+2)  (2n+l) 
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TABLE  II 

Integrals  in  Eg.  (30)  and  (31) 


Ju(ue-u)d£  =  J^2  +  J2$  +  J3 

^(ue-u)d£  =  J  4ip  +  j5 
|v(ue-U)d  =  jy  +  j7lfJ  +  j8 

^ve“v^d?  =  J18^  +  J19 

J(ve2-v2)dC  =  Jgip2  +  J1Q^  +  jn 

|v(ve-v)d5  =  J9»2  +  J12„  +  J13 

|u(ve-v)d?  =  J6ip2  +  J14^  +  J15 

J(ueve-uv)d?  =  j6^,2  +  J16ip  +  Jl7 
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TABLE  III 

Expressions  for  the  Jj 


J4 

J5 

J6 

J7 

J8 

J9 

J10 

J11 

Jl2 


J16 

J17 

J18 


=  A2G.^sin2a 

=  A  [-  (AG1~BG2)  cos  a  sin  a  -  2CG5sin2a] 

=  A(AG0+B) G^cos2a  +  C (-AG3+BG^) cos  a  sin  a  +  C2GgSin2a 
=  AG2sin  a 

=  AG7C0S  a  +  CG^sin  a 
=  -  A  G^  cos  a  sm  a 

=  A  [(AG^3+BG2) sin2a  +  AG34cos2a  +  2CG5cos  a  sin  a] 

=  [A(AG-^2+bg7)  “C2Gg  ]cos  a  sin  a  +  C (AGg+BG4) sin2a  -  ACG^Qcos2a 
=  A2G11cos2a 

=  A  [AG^- (qe+B) G2 ]cos  a  sin  a  -  2ACGgCOs2a 

=  A  [AG12+ (qe+B)  G?  ]sin2a  +  C  [AG3-(qe+B)G4]cos  a  sin  a  +  C2GgCOS2a 
=  -A  [(2AG-^3+2G2B-qeG2)  cos  a  sin  a  +  2  CGgCOS2a] 

=  [qe(2AG9+B)  -  (A2Gg+4G9AB+B2)  ]sin2a 

+C  [qgG^-2 (AGg+BG^ ) ]cos  a  sin  a  +  C2GgCOS2a 
=  A  [- (AG13+BG2) cos2a  +  (AG13+BG2-qeG2) sin2a  +  2CG5cos  a  sin  a] 

=  [qe(2AGg+B)  -  (A2Gg+4ABGg+B2) -C2Gg ]cos  a  sin  a 

— C (AGg+BG4 ) cos2a  +  C  (AGg+BG4~ q^G^) sin  a 
=  A  [- (AG^3+BG2)  (cos2a-sin2a)  +  2CG5sin  a  cos  a] 

=  [qe2- (A2Gg+4ABGg+B2) -C2Gg ]cos  a  sin  a  -  C (AGg+BG^) (cos2a-sin2a) 
=  -AG2cos  a 


Jjg  =  AGySin  a  -  CG4COS  a 


I 


r* 

? 


NOLTR  72-80 
TABLE  IV 

Program  for  Computation  of  Magnus  Force  and  Moment 

For  Half-Ellipsoid 


Column  I  are  symbols  that  occur  in  program 

Column  II  are  symbols  in  symbol  list  (page  vi)  that  corresi  r-3 
symbols  in  Column  1 


A 

AK 

AJN 

Alpha 

B 

C 

COS  A 

COSMA 

COSPH 

COST 

DELTA 

DELTAI 

DELTAZ 

DELTX 

DDDX 

DDSTAR 

DDSTDX 

DRDS 

DSTARI 

DTHDX 


II 


A 

k 

Ji 

a 

B 

C 


TT 


Cos  [  ct - ] 

180 

Cos  a 
cos  <p 

COS  0 
6 

Ss 

os  for  «  f  0,  u  =  0 
Ax 

dx 

A  6* 

d6*/dx 

dr 

_o 

ds 

6* 

s 

de 

dx 


T  IV- 1 


init 


DTHOLD 
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DUEXSP 

EN 


d0 

at  (x-&) 


s 


n 


FN 


G. 

1 


G 


K 


1 


GAMMA 


n 


H 


K 


2 


P 


P 


PDT 

PDX 

PJNT 

PJNX 

PPT 

PPX 

PRT 

PRX 

PRX1 


36 

9x 


dip 


3x 


u)r„ 


Vco 

3r 

o 

“3x 

1_  ffo 

ro  9x 


T  IV-2 
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51 

SII 

UET 

•UEX 

5  VET 


QDQN 

QDQS 

R 

REL 

RHO 

RPDSTR 

RUC 

SIN  A 

SMA 

SNSMA 

SINT 


% 


3u 


e 

3x 


3V, 


39 

3V 


e 

3x 


+  V 


„  3qe 
qe  ~i 

qe 

3S 


Reynolds  number 


VooL 


r  +  5* 

G 


r  u  cos  d> 
o  e 


sin[aI$0] 


a 

Sin  a 


Sin  0 


T  IV- 3 


T 

THETA 

TWQE 

TWRT 

TWXR 

URATIO 

V  RATIO 

XEN 

XP 

XSTAG 
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t 

0 


Value  of  n  for  x=x 

s 

Values  of  x  chosen  for  printout 
before  run  on  computer 
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